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Abstract 



We complete the construction of the renormalization group equation (RGE) 
for the Color Glass Condenstate begun in Paper I. This is the equation which 
governs the evolution with rapidity of the statistical weight function for the 
color glass field. The coefficients in this equation — one- loop real and virtual 
contributions — are computed explicitly, to all orders in the color glass field. 
The resulting RGE can be interpreted as the imaginary-time evolution equa- 
tion, with rapidity as the "imaginary time", for a quantum field theory in 
two spatial dimensions. In the weak field limit it reduces to the BFKL equa- 
tion. In the general non-linear case, it is equivalent to an equation by Weigert 
which summarizes in functional form the evolution equations for Wilson line 
operators previously derived by Balitsky and Kovchegov. 



1 



I. INTRODUCTION 



Hadronic scattering at high energy, or small Bjorken's x, uncovers a novel regime 
of QCD where the coupling is small (a s <C 1) but the parton densities are so large that 
conventional perturbation theory breaks down, via strong non-linear effects 0-lO|. In a 



previous paper [PH| , to be referred to as "Paper I" throughout this text, we have outlined 
the construction of an effective theory which is well suited to study the non-linear phenomena 
at small x and has a transparent physical interpretation: It portrays the gluon component of 
the hadron wavefunction (the relevant component at small x) as a Color Glass Condensate 
(CGC). This is a multiparticle quantum state with high occupation numbers, but to the 
accuracy of interest it can be represented as a classical stochastic color field with a probability 
law determined by a functional Fokker-Planck equation. 

The latter is a renormalization group equation (RGE) which shows how to construct 
the effective theory by integrating out quantum fluctuations in the background of a strong 
color field (the CGC). The non-linear effects included in the RGE via this background field 
describe interactions among the gluons produced in the quantum evolution towards small x. 
This leads to a non-linear generalization of the BFKL equation |T2|] whose general structure 



has been originally identified in Refs. [|13 , in an effort to give a field-theoretical justification 



to the McLerran-Venugopalan model for the gluon distribution of a large nucleus 0,n-TB 
But previous attempts |17| to compute the coefficients in this equation (beyond the linear, 
or BFKL, approximation |T3|]) have suffered from technical complications and, moreover, 
appear to be inconsistent with the results obtained from other perturbative approaches 



In this paper, we continue the analysis in Paper I and compute explicitly the coefficients 
in the RGE alluded to above. As recently reported in Ref. [pT| , the non-linear effects encoded 



in these coefficients lead to gluon saturation, that is, to a limitation on the maximum gluon 
density in the hadron wavefunction as x — > 0. In contrast to the linear BFKL equation, 
which predicts the exponential growth of the gluon distribution function with r = ln(l/x), 
our corresponding result grows only linearly, and thus respects the unitarity bounds. This is 
consistent with a previous result by Mueller J7j] , and also with some recent analyses [F22J-f2^fT^] 
of the Balitsky-Kovchegov equations |18|.19[|, which are encoded too in our RGE |25| (see 



also Sect. |VB| below). 

This paper is quite technical and relies heavily on the results of Paper I. To facilitate 
its reading and prepare the calculations to follow, it is convenient to summarize, in the 
remaining part of this Introduction, the general structure of the effective theory and its 
quantum evolution. In this introductory discussion, we shall follow the presentation in Paper 
I, to which the reader may refer for more details, but we shall also anticipate some of the 
results to be obtained later in this paper. The remaining sections are organized as follows: In 
Section II we present the explicit calculation of the "real correction', i.e., the charge-charge 
correlator x which plays the role of the diffusion kernel in the functional RGE. In Section III 
we do the same for the "virtual correction", i.e., the one-point function a which plays the 
role of a force term. In Section IV, we present the final result for the RGE, and discuss its 
general structure and some of its remarkable properties. We emphasize, in particular, the 
Hamiltonian structure of this equation, and its relation with a similar equation by Weigert 
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2Cfl . In Section V, we derive evolution equations for observables from the general RGE. In the 
weak field limit, we thus recover the BFKL equation for the unintegrated gluon distribution. 
In the general non-linear case, we obtain the coupled evolution equations for Wilson-line 
correlators originally derived by Balitsky and (in the large N c limit) also by Kovchegov. In 
Section VI we present conclusions and perspectives. Some explicit calculations, as well as the 
presentation of the background field gluon propagator, are deferred to the four Appendices. 



A. The effective theory for the CGC 



The effective theory applies to gluon correlations in the hadron wavefunction as mea- 
sured in deep inelastic scattering at small Bjorken's x. It is formulated in the hadron infinite 
momentum frame, where small x corresponds to soft longitudinal momenta^ k + = xP + , with 
P + the hadron momentum, and i < 1. The main physical assumption, which is a posteriori 
verified in the construction of the effective theory, is that the "fast" partons in the hadron 
wavefunction (i.e., the excitations with relatively large longitudinal momenta p + ^> k + ) can 
be replaced, as far as their effects on the soft correlation functions are concerned, by a clas- 
sical random color source p a (x), whose gross properties are determined by the kinematics. 

The separation of scales in longitudinal momenta (p + ^> k + ) implies a corresponding 
separation in (light-cone) energies (p~ ~ p\/2p + <C k~ ~ kj_/2k + ), and therefore also in 
(light-cone) time scales: The lifetime Ax + ~ l/k~ cx k + of the soft gluons is much shorter 
than the characteristic time scale ~ oc p + for the dynamics of the fast partons. Thus, 
the latter appear to the soft gluons as nearly on-shell colored particles which propagate 
on the light-cone (z ~ t, or x~ ~ 0) with large p + momenta. The associated color charge 
density p a (x) is therefore static (i.e., independent of x + ), localized near the light-cone (within 
a small distance Ax - ~ "C l/k + ), and random (since this is the instantaneous color 
charge in the hadron "seen" by the soft gluons at the arbitrary time of their emission). 
The correlations of p are encoded in a gauge- invariant weight function W T [p}. This is the 
probability density for having a color charge distribution with density p a {x~ , XjJ, normalized 
as: 

fv P W T \p] = l. (1.1) 

Note that we use the momentum-space rapidity r = ln(P + /k + ) = ln(l/a;) to indicate the 
dependence of the weight function upon the soft scale k + . 

Thus, in this effective theory, the equal-time[| gluon correlation functions at the scale 
k + = xP + are obtained as (with x = (x _ ,xj_)): 



throughout, we use light-cone vector notations, that is, = (v + ,v~,v±), with v + = 
(l/v / 2)(^° + v 3 ), v~ = (l/y/2)(v° - v 3 ), and vj_ = (v 1 ^ 2 ). The dot product reads k ■ x = 
k~x + + k + x~ — kj_ • xj_. The hadron four-momentum reads = (-P + , 0, 0j_), with P + — > oo. 

2 Only equal-time correlators are needed, since these are the only ones to be measured by a 
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(Al(x + ,x)Al(x + ,y)---) T = Jvp W T [p}A l a (S)Ai(y)---, (1.2) 

where A a = .44 [p] is the solution to the classical Yang-Mills equations with source p a : 

(D v F"%(x) = c^ + p a (x), (1.3) 

in the light-cone (LC) gauge = 0, which is the gauge which allows for the most direct 
contact with the gauge-invariant physical quantities P,|TT||. For instance, the gluon distribu- 
tion function (= the total number of gluons per unit rapidity with longitudinal momentum 
k + = xP + and transverse momentum k\ < Q 2 ) is obtained as P,|TT 



xG(x,Q 2 ) = U^e(Q 2 -kl)(\^(k)\ 2 ) . (1.4) 



7T J (2?r) 

where = d + A a is the electric field associated to the classical solution .44 [p], and k = 
(fc+, kj.) with k+ = xP+ = P+e~ T . 

To deduce explicit expressions for these classical fields, it is preferable to express the 
LC-gauge solution A a in terms of color source p a in the covariant gauge d^A^ = (COV- 
gauge). This is possible since both the measure and the weight function in the functional 
integral (|1.2j ) are gauge-invariant, so that the classical average can be done equally well by 
integrating over the COV-gauge p. In terms of this latter, the classical solution A%[p] is 
known explicitly and reads: = (the gauge condition), A~ = 0, and (in matrix 



notations appropriate for the adjoint representation: p = p a T a , etc.) 



A i (x) = -U(x)d i UUx), (1.5a) 
9 

W(x~, x±) = P exp J dz~ a(z~~ , , (1.5b) 

-Vj.a(x) = p(x). (1.5c) 

This is the gauge-transform, with gauge function U(x), of the corresponding solution in the 
COV-gauge, which has only one non-trivial component: A% = 5^ + a a {x). All the fields above 
are static, i.e., independent of x + . Since p, and therefore a, are localized at Ax~ < l/k + , 
the associated field A appears effectively as a ^-function : 

A*(x-,x ± ) w 0(x-)-V(&V*) = 0(x-)A*(x ± ), (1.6) 
9 

to any probe with a much lower longitudinal resolution (i.e., with momenta q + k + ). On 
the same resolution scale: 



small-x external probe, which is absorbed almost instantaneously by the hadron. In eq. ( |1.2| ), 
it is understood that the fields A l a (x + , x) involve only Fourier modes with longitudinal momenta 
k+ = X P + . 
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rf(x-,x ± ) « 6{x-)V*{x±) + 9{-x-), F+\x) « 5(aT)^(x ± ). (1.7) 

In the equations above, V and are the asymptotic values of the respective gauge rotations 
as x" — ► oo : 

V^xj.) = Pexp jig y dzT a(#~, xjj | . (1.8) 

In practice, £7(x~,xj_) = V^(xj_) for any x~ ^> l/k + . 

Note that the solution to the classical non-linear equation is known exactly, which 
makes this approach particularly convenient to study the non-linear physics at small x. The 
classical approximation should be indeed appropriate in this regime, which is characterized 
by weak coupling and large occupation numbers. On the other hand, the weight function 
W r [p] is obtained via a quantum calculation in which the quantum fluctuations with mo- 
menta k + ^ p + ^ P + are integrated out in layers of p + , and in the background of the 
classical fields A 1 [p~5| , p~3] , [T/T]| . This calculation captures the basic mechanism leading to large 
gluon densities — namely, the BFKL-type of evolution towards small x — , while also in- 
cluding non-linear effects (i.e., rescatterings among the produced gluons) via the background 
fields. In the saturation regime, where one expects color fields as strong as A 1 ~ 1/g (cor- 
responding, via eq. (|1.4|), to gluon densities of order l/a s HJT5||7|,fH|1), the background field 
calculation must be carried out exactly, i.e., to all orders in A 1 . 



B. The quantum evolution of the effective theory 



To describe the quantum evolution, it is convenient to introduce an arbitrary separa- 
tion scale A + such as k + ^ A + <C P + , and assume that the "fast" quantum modes with 
momenta p + ^> A + have been already integrated out and replaced, to the accuracy of in- 
terest, by a classical color source p a (x) with weight function Wa[p]- On the other hand, 
the "soft" gluons with momenta p + < A + are still explicitly present in the theory. The 
correlation functions at the soft scale k + are then obtained from the following generating 
functional : 

r { f A VA^S(A+) e iS[A,p]-ifj-A] 

Z[J] = I v ? WM { s ^ml t) ^ ) ■ (L9) 

where the "external current" J£ is just a device to generate Green's functions via functional 
differentiations, and should not be confused with the physical source p. Eq. (|1.9|) is written 
fully in the LC gauge (in particular, p is the LC-gauge color source), and involves two func- 
tional integrals: (a) a quantum path integral over the soft gluon fields A 11 , which generates 
(p-dependent) quantum expectation values at fixed p, e.g., 

C™'(*M'to»H = (1.10) 
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and (b) a classical average over p. As discussed in Paper I, eq. ( |1.9|) has the typical structure 
to describe correlations in a glass. 

The action S[A, p] in eqs. ( |1.9| ) and ( |1.10| ) reads as follows [13|] : 
S[A,p] = - J d A x- A F; v F^ + -^Jd 3 xTr{p(x)W[A-}(x)} = S YM + S w , (1.11) 

where VTfA"] is a Wilson line in the temporal direction: 



W[A-](x) = P exp 



ig / dx + A (x) 



1.12) 



In the classical, or saddle point, approximation 8S/8A^ = 0, the action ( |1.11| ) generates 
the desired equations of motion, that is, eqs. (|L3l ) with A~ = 0. This shows that the 
classical field in eqs. ( |1.5| ) is the tree-level background field in which propagate the quantum 
fluctuations. The general non-linear structure of Sw in eq. ( p.,12| ), which plays a role only for 
the quantum corrections (in that it generates new interaction vertices; cf. Sect. [TT] below), 
reflects our eikonal approximation for the interaction between fast particles moving in the 
plus direction (here represented by p) and the comparatively slow gluon fields. 

In the above formulae, the intermediate scale A + enters at two levels: as an upper 
cutoff on the longitudinal momenta p + of the quantum gluons, and in the weight function 
Wa[p] for the classical source. Of course, the final results for correlation functions at the 
scale k + must be independent of this arbitrary scale : 



A" 



dZ[J] 
dA+ 



0. 



;i.i3) 



This constraint can be formulated as a renormalization group equation (RGE) for Wa[p] 
which governs its evolution with decreasing A + . 

The initial conditions for this evolution are determined by the properties of the 
hadronic matter at A + ~ P + . These are not really under control in perturbation the- 
ory, but one can try to rely on some non-perturbative model, like the valence quark model. 
The initial conditions might be under better control in the high density environment of a 
very large nucleus, since we expect the coupling constant to be weaker at high density. The 
indeterminancy of evolution equations due to initial conditions is hardly a new problem in 
QCD, since the DGLAP and BFKL equations are both limited by such uncertainty. In 
these cases, and we also hope hereQ, the solution of the evolution equation for arbitrarily 
high energies is universal and its generic properties are largerly independent of the initial 
conditions. 

Starting with these initial conditions, one then proceeds with the quantum evolution 
down to the soft scale of interest A + ~ xP + <C P + . In this process, the original source 



3 The approximate solutions to the RGE recently found in Ref. [21] appear to confirm this 
expectation. 
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at the scale P + gets dressed by the quantum fluctuations with momenta A + < \p + \ < P + . 
We treat this process in perturbation theory, in the "leading logarithmic approximation" 
(LLA)[j — i.e., by retaining only the terms enhanced by the large logarithm ln(l/x), to all 
orders in (a s ln(l/x)) n ) — , which, by itself, is also the accuracy of the BFKL equation, 
but we go beyond BFKL in that we resum also finite density effects, which are expected to 
become increasingly important as we go to smaller and smaller A + (or Bjorken's x). We do 
that by performing a background field calculation, that is, by integrating out the quantum 
fluctuations at one step in the background of the classical field generated by the color source 
at the previous step, with the background field simulating (via its correlations) the finite 
density effects in the system. 

To be more explicit, let us describe one step in this renormalization group procedure 
in some detail. Assume that we know the effective theory at the scale A + - - as specified 
by the corresponding weight function Wa[p] = Wr[p], with r = ln(P + /A + ) — , and we are 
interested in correlations at the softer scale k + ~ bA + with b <C 1, but such as a s ln(l/6) < 1, 
for perturbation theory in powers of a s ln(l/b) to make sense. Our purpose is to construct 
the new weight function W&a[/o] = Wt+At[p], with At = ln(l/6), which would determine 
the gluon correlations at this softer scale. As compared to W r [/>], this new weight function 
must include also the quantum effects induced by the "semi-fast" gluons with longitudinal 
momenta in the strip 

bA + < \p + \ < A+. (1.14) 

To compute these effects, it is convenient to decompose the gluon field in eq. ( |1.9| ) as follows 

A* = A*\p] + a* + 5A*, (1.15) 



where A% is the tree-level field, represents the semi-fast fluctuations, and 5A^ refers 
to the remaining modes with \p + \ < bA + . By integrating out the fields a^ 4 , some new 
correlations are induced at the soft scale 6A + via the coupling 5A~5p c , where 5p c is the 
color charge of the semi-fast gluons. These correlations have to be computed to lowest order 
in a s ln(l/6), but to all orders in the background fields A % . This is essentially an one-loop 
calculation, but with the exact background field propagator of the semi-fast gluons. The 
new correlations are eventually absorbed into the functional change AW = W t +At — W T in 
the weight function. Since AW oc At, this evolution is most conveniently formulated as a 
(functional) renormalization group equation for W T [p], which reads ]13| , [ll] 

= S \\ x , * d W rX ? y ] - jr^—dW^]) , (1.16) 
Ot s \2 5p-(x ± )6p b T (y ± ) 1 TAxyS 5p a T (x ± ) [ lJ J 

in compact notation where repeated color indices (and coordinates) are understood to be 
summed (integrated) over. The coefficients a" = cr a (xj_) and Xxy = Xab(%±,V±) in the above 



4 Indeed, it is only to this accuracy that the assumed separation of scales in the problem is 
maintained by quantum corrections. 



equation are related to the 1-point and 2-point functions of the color charge Sp a (x) of the 
semi-fast gluons via the following relations: 



a s \n-a a (x±) 



dx (Sp a (x)) , 

a s \n^-Xab{x±,y±) = J dx" J dy~ (5p a (x + ,x) 5p b (x + ,y}) 



;i.i7) 



where the brackets denote the average over quantum fluctuations in the background of the 
tree-level color fields A 1 , as shown in eq. ( 1.10Q . 



In writing eq. ( |1.16| ), we have also anticipated the longitudinal structure of the quan- 
tum evolution, which will become manifest only after performing the quantum calculations 
in the next sections. Specifically, we shall see that the induced source (Sp a (x)) (= the 



correction to p generated by the gluons with p + in the strip ( |1.14j )) has support at 

1/A+ < x- < 1/(6A + ). 



;i.i8) 



By induction, we shall deduce that p a (x) (= the color source generated by the quantum 
evolution down to A + = P + e~ T ) has support atf] < x~ ^ x~ with 



x. 



1/A H 



x e 



1/P H 



1.19) 



Thus, the color source is built in layers of x~, with a one-to-one correspondence between the 
^"-coordinate of a given layer and the p + -momenta of the modes that have been integrated 
to generate that layer. When some new quantum modes, with rapidities r < r' < r+Ar, are 
integrated out, the preexisting color source at < x~~ < x~ remains unchanged, but some 
new source is added in the interval ( |1.18|) . Because of that, AW = W t +At — W T involves 
only the change in p a within that last interval. In the continuum limit, this generates the 
functional derivatives of W T with respect to p a (x) at x~ = x~ only, that is, the derivatives 
with respect to 



p a (x~ 



X„ 



x±), 



;i.2o) 



as shown in eq. ( |1.16| ). 

On the other hand, given the separation of scales in the problem, the detailed longi- 
tudinal structure of the quantum effects should not be too important, and this too has been 
anticipated in writing eqs. ( 1.16 )— ( 1.17 ): In the same way as the original source with support 
at x" ^ 1/A + appears effectively as a 5-function to the semi-fast gluons (with wavelenghts 
Ax" ~ l/p + 3> 1/A + ), the induced source (6p), although relatively delocalized as com- 
pared to p (cf. eq. ( |1.18| )), appears to the soft gluons with k + ^ 6A + as a rather sharp color 
distribution around x~ ~ x~. Thus, the soft gluons can probe only the two-dimensional 



5 The fact that the source has support at positive x , rather than at generic x with \x 



< 



will be seen to be related to our specific gauge-fixing prescription; cf. Sect. Ill A 



S 



color charge distribution in the transverse plane, as obtained after integrating out x , cf. 
eq. Ql.17 ). This explains the 2- dimensional structure of the coefficients a and x i n the RGE 



(El)- 

Eq. (|1.16|) has the structure of the diffusion equation: It is a second-order (functional) 



differential equation whose r.h.s. is a total derivative (as necessary to conserve the total 
probability; cf. Sect. |I V A|) . Thus, this equation describes the quantum evolution as a 



diffusion (with diffusion "time" r) of the probability density W T [p] in the functional space 
spanned by p a (x~,x±). For this equation to be useful, its coefficients a and x must be 
known explicitly as functionals of p. It is therefore more convenient to use the COV-gauge 
source p a , or the associated Coulomb field a a , as the functional variable to be averaged 
over. Indeed, a and % depend upon the color source via the classical field A\ which is most 
directly related to a a , cf. eqs. (|1.5|). 



As explained in Ref. [11 1, when going from the (background) LC-gauge to the COV- 
gauge, the induced color charge a acquires a new contribution, in addition to the color 
rotation with matrix W. This is so because the transformation between the two gauges 
depends upon the charge content in the problem: this was p at the initial scale A + , thus 
giving a gauge transformation £/^[p], but it becomes p + 5p, with fluctuating 5p, at the 
new scale 6A + , thus inducing a fluctuating component in the corresponding gauge function 
U^lp + Sp]. After averaging out the quantum fluctuations, one is left with a RGE for W T [p] 
which is formally similar to eq. (|1.16| ), but with p — * p, x ~ > X, an d cr — ► a, where: 

Xab(x±,y±) = V^ c (x±) Xcd(x±, y±) V db (y±), (1.21a) 

&a(x±) = KU X -L) cr b(x±) ~ 8<7 a (x±), (1.21b) 

Sa a (x ± ) = |r 6c / d 2 y ± Xcb(x ± ,y ± ) (y ± \ \x ± ) , (1.21c) 

with V and as defined in eq. ( |1.8|) . The correction —5a, to be referred to as the "classical 
polarization" , is the result of the quantum evolution of the gauge transformation itself. 

It turns out that the RGE is most conveniently written as an equation for W r [a] = 
W T [p = — Vj_a], in which case it reads: 

^ = a - { ^(S^y,) lW ^ - ^k) lWr "' ] } • (L22) 

where a®(x±) = a a (x~ = x~,x±) (cf. eq. ( |1.20[ )), and 

v a (x ± ) = J d 2 z ± (x ± \ \z ± ) d a ( Zl _) } (1.23a) 

V ab (x ± ,y ± )= Jd 2 z ± fd\ ± (x ± \-^\z ± )x ab (z ± ,u ± )(u ± \-^r\y ± ). (1.23b) 

This form has, in particular, the advantage to eliminate any explicit reference to the arbitrary 
infrared cutoff p, which enters the solution to eq. (|1 .5cj ) for a : 



d 2 

a{x) = - I ln(|x ± -y±\p) p(x~,y±). (1.24) 
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It is our ultimate goal in this paper to compute explicitly the quantities r] and v as functionals 
of the field a, and thus completely specify the RGE ( |1.22|) . 



C. Feynman rules for \ and a 

The basic quantum calculation is that of the quantities x an d o in the LC gauge. 
The necessary Feynman rules can be readily derived from the action (|1 . 1 1| ) |y|[lT], and are 



summarized here for later reference. Eqs. ( |1.17| ) involve the color charge Sp a of the semi-fast 
fields a^. To the order of interest, this reads 

5p a (x)=5p^(x) + 5pf\x), (1.25) 

where Sp^ is linear in the fluctuations a^, while 8p^ is quadratic^: 

8pM(x) = -2igF+ i (x)a ic (x) + 

+gp ac (x) J dy + (x + \PV^ z \y + )a c -(y + ,x), (1.26) 
5p^\x)=gr bc [d + a\{xM{x) 

- (g 2 /N c )p\x) J dy + a~ c {y + ,x) J dz + a~ d {z + ,x) 

x \e(z + - y + )6(y + - x + )Tr (T a T c T d T b ) 
+ 6{x + - z + )6{z + -y + )Ti (T a T h T c T d ) 

+ 6(z + - x + )6{x + - y + )Ti (T a T d T fe T c )|. (1.27) 

In the right hand sides of these equations, the terms involving a l c come from the three-gluon 
vertex in Sym, while the terms involving a~ come from the two- and three-point vertices in 
Sw- It is understood here that only the soft modes with k + ^ bA + are to be kept in the 
products of fields. 

Consider first the following charge-charge correlatorf] which enters eq. ( |1.17|) for \ '■ 



Xab(x,y) = (5p a (x + } x) 5p b (x + ,y)) . (1.28) 

By time homogeneity (recall that the background fields A 1 are static), this equal-time 2-point 
function is actually independent of time. When evaluating \ within our present accuracy, it 



6 To compare with the corresponding formulae in Paper I (see, e.g., eqs. (4.41)^(4.42) in Ref. 
|i"T|l ), please note a change in our normalization conventions: Here, the various fields preserve their 
natural normalization fixed by the action ( |l.ll| ). By contrast, in Paper I we have rescaled the 
classical fields and sources by a factor 1/g. 

7 Our notations are such that \ an d X or o an d o differ just by a factor of a s ln(l/6). 
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is sufficient to keep only the terms of first order in fluctuations: 5p ~ 5p^\ These are given 
by eq. ( |1.26[ ) which implies (in condensed notations, where the PV prescription in l/p~ and 



the condition x + = y + are implicit): 



1 \ a /_. 1 x bs 



X ab {x,y) = g 2 ((-2iF +l a l + p—a-j [2ia % T +l + a' — p) ). (1.29) 



This involves the propagator 

iG^(x,y)[A] = (Ta^x)a u b (y)) (1.30) 

of the semi-fast gluons in the background of the tree-level fields A\ and in the LC-gauge 
(a + = 0). Given the specific structure of the background field ( |1.5| ), it has been possible to 
construct this propagator exactly, i.e., to all orders in A 1 |39|j40|JTT[| . The resulting expression, 
to be extensively used in what follows, is presented in Appendix A. 

The construction of the propagator in Refs. [^,||y.|ll|] has relied in an essential way 
on the separation of scales in the problem: Because of their low p + momenta (p + <C A + ), 
the semi-fast gluons a M are unable to discriminate the longitudinal structure of the source, 
but rather "see" this as a ^-function at x~ = 0. One is thus reduced to the problem of 
the scattering off a £-type potential, whose solution is known exactly. But this also means 
that, strictly speaking, the propagator ( |1.30| ) is known only far away from the support of 
the source (for \x~\ ^> 1/A + ), where A" takes the approximate form in eq. (|1.6| ). It is thus 
an important self-consistency check on the calculations to follow to verify that \ and a are 
indeed insensitive to the internal structure of the source. 

Still for the construction of the propagator, it has been more convenient to impose the 
strip restriction (1T4 ) on the LC-energy p~ , rather than on the longitudinal momentum p + 



||1 1|| . Indeed, to LLA, the quantum effects are due to nearly on-shell gluons, for which the 
constraint fll.l4| ) is equivalent to the following constraint on p~ : 

A" < \p-\ < A"/6, (1.31) 

where A~ = Qj_/2A + , and Q± is some typical transverse momentum. (Recall that we 
assume all the transverse momenta to be comparable; thus, to LLA, it makes no difference 
what is the precise value of Q±.) But the condition ( |1.31| ) is better adapted to the present 
background field problem, where we have homogeneity in time (so that we can work in the 
^-representation), but strong inhomogeneity in x~ (cf. eqs. ( |1.6| ) and (|L~7|)). 

A final subtlety refers to the choice of a gauge condition in the LC-gauge propagator: 
Recall that, even after imposing the condition A + = 0, one has still the freedom to perform 
^'-independent gauge transformations. This ambiguity shows up as an unphysical pole at 
p + = in the gluon propagator. Since in our calculations with strip restriction on p~, p + is 
not restricted anylonger, an ie prescription is needed to regulate this pole. This is important 
since, as our calculations show, the final results for \ and a are generally dependent upon 
this ie prescription. In fact, the discrepancy between our results below and those reported 
in Refs. ]17| may be partially attributed to using different gauge-fixing prescriptions. 
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For consistency with the retarded boundary conditions (A 1 — > for x~ — > — oo) im- 
posed on the classical solution (|1.5| ), we shall use a retarded prescription also in the 
LC-gauge propagator of the quantum fluctuations (see Appendix A). With this choice, the 
induced source (5p) appears to have support only at positive x~, as shown eq. ( |1.18|) . We 



have also verified that, by using advanced prescriptions (both in the classical and the quan- 
tum calculations), the support of (Sp) would be rather at negative x~ (but with \x~\ still 
constrained as in eq. (|1.18 )). In both cases, however, the same results are finally obtained 



for x and a after integrating out the (prescription-dependent) longitudinal structure. Thus, 
the RGE comes up the same with both retarded and advanced prescriptions. On the other 
hand, we have not been able to give a sense to calculations with other prescriptions, like 
principal value or Leibbrandt-Mandelstam. (For instance, when computed with a principal 
value prescription, the coefficients a and x appear to be sensitive to the internal structure 
of the source in x~ , which contradicts the assumed separation of scales.) 

We conclude this introductory section by giving explicit expressions for x an d cr in 
terms of the background field propagator G^(x,y) of the semi-fast gluons. For x, we use 
eqs. ( |1.17| ) and ( |1.29| ) to deduce that 



a s ln^Xafe(x_L,y_L) = J dx J dy Xab(x,y), (1.32) 

with (in matrix notations) 

- 2 X (x , y) = i 2Tf (x\G*\y) VF* + 2J*< (xp' -L|y) p y 

- Px ( x \^-G' i \y)2T^ i + ip^xl-^-G-^^py. (1.33) 

The equal-time limit is implicit here; it is achieved by taking y + = x + + e within the 
time-ordered propagator (|1.30| ) [ |TT[] . Diagramatically, all the above contributions to x are 
represented by tree-like Feynman graphs (no loops), with vertices proportional to p or jF +i 
(see Fig. |]). Since, physically, these are quantum corrections associated with the emission 
of a real (semi-fast) gluon, we shall sometimes refer to x or x as the real correction. (In the 
weak field approximation, x is responsible for the real piece of the BFKL kernel p!3|JTT| ; see 



also Sect. VA below. 
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FIGURES 




a) 






C) 



d) 



FIG. 1. Feynman diagrams for the four contributions to x given in eq. ( 1.33| ), 



To the same accuracy, the induced source (5p) involves only the terms in 5p of second 
order in the fluctuations : 



cr. 



a (x) = (Sp a (x)) = (5tf\x)), 



with 5p^ given by eq. (|1.27Q . This yields 

a s lnia a (x_L) = J dx~ Tr (T a a(x)), 

where: 



1.34) 



1.35) 



a(x) = -gd+G u (x, y) + ig 2 p(x)(x 

x=y \ 



1 G 1 



id~ id~ 



a Ax) + a 2 {x). 



1.36) 



In writing &2 as above, we have used compact but formal notations for the second contri- 
bution to 5p^ in eq. ( |1.27| ), which is non-local in time. Diagramatically, the two terms in 
eq. ( |1.36D are represented by the one-loop diagrams displayed in Fig. 0. These are vertex 
and self-energy corrections which in the weak field limit (i.e., to linear order in p) generate 
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the virtual piece of the BFKL kernel [|T^,|TT|]. Accordingly, we shall refer to a or a as the 

virtual correction. 





a) b) c) d) 

FIG. 2. Feynman diagrams for d\ (a) and &i (b,c,d). The wavy line with a blob denotes the 
background field propagator of the semi-fast gluons; the continuous line represents the source p; 



the precise vertices can be read off eq. (1.27). 



II. THE INDUCED CHARGE-CHARGE CORRELATOR, OR "REAL 

CORRECTION" 



We are now prepared for the explicit calculation of the coefficients in the RGE, to 
be completed in this and the following section. The one-loop quantum calculation will be 
performed fully in the LC-gauge (a+ = = 0), by using the retarded ie prescription 
discussed in Sect. [I C| . The corresponding background field propagator is given explicitly in 
Appendix A. The ensuing expressions for \ an d cr, which are gauge- covariant functionals of 
the background fields, will be then rotated to the COV-gauge for the background fields (cf. 
eq. (|1.21 )), which is the only gauge which allows for explicit non-linear calculations. 



In particular, the calculation of the "real correction" in this section will be organized 
as follows: We shall first compute x(x, y) by evaluating the matrix elements in the r.h.s. of 
eq. ( |1.33| ), then we shall derive x{ x ±iV±) by integrating out the longitudinal structure of 
x(x, y), cf. eq. ( |1.32| ). The result will be subsequently rotated to the (background) COV- 
gauge, to give x(xj_, y±), cf. eq. ( |1.21| ). Finally, the coefficient r/(x±, y±) in the RGE for 
PF r [a] will be obtained according to eq. ( |1.23b| ). 



A. Non-linear effects and their dependence upon the gauge-fixing prescription 



By substituting the expressions (|A12| ) for the gluon propagator in eq. ( |1.33j ) for x{%, y), 
one obtains, after simple algebra, 

^ X (x,y)=i2^ i G i \x,y) 
9 
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= T^i^iT) + Xa(^, j/) + Xs{x,y)), (2.1) 

where the equal-time limit y + = x + + e is implicit, and Xi, X2, and x 3 refer, respectively, to 
the terms in the first, second and third line. 

In these expressions, T> % = d % — igA 1 and = + igA^ (with the derivative 
acting on the function on its left) are covariant derivatives constructed with the background 
field A 1 . Furthermore, G^(x,y) is the gluon propagator in the temporal gauge d~ = 0, and 
is presented in Appendix A. 3. It is a non-linear functional of the background field a(x), via 
the Wilson lines V and (see, e.g., eqs. (|A8|)- (|A~9|) ). What is however remarkable, and 



will be demonstrated in what follows, is that all the non-linear effects encoded in G^ v drop 
out in the calculation of the matrix elements in eq. ( [2.1|) . That is, the final result for x(x, y) 
is the same as obtained by evaluating the r.h.s. of eq. (|2.1| ) with the free temporal-gauge 
propagator Gq U . This simplification is a consequence of our specific ie prescription in the 
LC-gauge propagator, as we explain now. Consider the following matrix element: 

(2-2) 

which enters the r.h.s. of eq. (|2.1| ). Our "retarded" prescription in G^ u implies that, on the 
left of G ++ in eq. ( |2.2| ), = id + +ie should be retarded, while on the right -r^+ = id +_ ie 
should rather be advanced (cf. eqs. ( |A13| )- (|A14|) ). Thus, eq. (|2.2|) is the same as 



dz l I dz 2 (x I 1 , \z x ) (x + ,z 1 ,x±\ G ++ \y + ,z 2 ,y±) (z 2 \ 1 , \y ) 
io + + ie io + — it 

dzi / dz 2 6(x~ - zf ) G ++ (z^,z 2 ) 9(y~ - z 2 ) , (2.3) 



where in the second line only the longitudinal coordinates have been shown explicitly. In 
eq. (|2.1| ), this matrix element is sandwiched between p{x) and p{y). Since p(x) is localized at 
smallx - (x _ ^ x~, cf. eq. (|1.19|) ), while G tJ,v (x, y) is relatively slowly varying as a function 
of x~ and y~ (since this is the propagator of the semi-fast gluons, with p + A + ), one can 
effectively replace x~ ~ and y~ ~ in eq. (|2.3| ) [and everywhere else in eq. (|2.1| ); recall 
that the electric field T +l is as localized as p\. This gives, for x + = y + , 

(x~~0, Xl _\^G ++ ^- \y~ ~0, y ± ) = J dz^ J dz 2 9(-z^)G ++ (z^ z 2 )9(-z 2 ) , (2.4) 

which shows that both external points , z 2 in G are negative, so G must be non-crossing. 
(Recall that "crossing" and "non-crossing" refer to whether the gluon has propagated or not 
accross the surface at x~ = 0, where the color source is located; see Sect. 6 of Paper I and 
Appendix A in this paper.) Moreover, the non-crossing piece of G(zi , z 2 ) at negative z~\ , z 2 
is the same as the corresponding piece of the free propagator. Thus, the matrix element 
( p.4[) singles out that particular piece of G which is not affected by the background field. 
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A similar conclusion applies to the other terms in the r.h.s. of eq. (|2.1| ), which involve 
the following matrix elements 

(x\&\y), ( x \-L-G + *\y), ( x \G + -^\y), (2.5) 

evaluated at x~ ~ y~ ~ and y + = x + + e. For instance, in the first term above, one can 
use the continuity of at x~ = and y~ = (cf. Appendix A) to approach these points 
from x~ < and y~ < 0, i.e., from the domain where G^ coincides with the corresponding 
free propagator Gq . We thus deduce that 

(0,x ± \G ij \0,y ± ) = (0,x ± \G%\0,y ± ), (2.6) 



and similarly for the other matrix elements in eq. 

For the previous arguments, it has been essential that the retarded prescription has 
been used systematically in the LC-gauge, both in the classical solution and in the quantum 
propagator: (a) The retarded boundary condition on the classical solution has insured that 
the background field A 1 is non- vanishing only at x~ > 0, or z < t. That is, the classical field 
sits behind its source, the (fast degrees of freedom of the) hadron, which is located at z = t. 
(b) The retarded ie prescription in the gluon propagator has implied that the semi-fast gluon 
exchanged within \ (cf. Fig. |]) is confined at negative x~, where there is no background 
field. Thus, this quantum gluon propagates freely from y to x, with x~ ~ y~ ~ 0. 

Note, however, that this property of a free propagation holds only for the temporal 
gauge gluon d M (with a" =0), which does not couple directly to the singular color source 
p [[LI]]. Because of that, its propagator G^ u {x,y) is continuous at x~ = and y~ = 0, 



and the associated non-linear effects drop out in the calculation of x, as argued before. By 
contrast, the light-cone gauge propagator G^(x,y) is sensitive also to the discontinuous 
gauge rotations at the end points x~ and y~~, via the covariant derivatives T> x and T>y 
(which technically enter via the gauge rotation from the temporal gauge to the LC gauge; 
cf. eq. ( |A11| )). As obvious on eq. (|2~T|), the non-linear effects associated with the fields 
within T> l x and T>y do subsist in the final result for \. Thus, while it is correct to replace 
the temporal gauge propagator G^ v in eq. (|2.1|) by its free counterpart, such a replacement 
would not be legitimate for the light-cone gauge propagator G^ u in eq. ( |1.33| ). 

By inspection of the previous arguments, it should be also clear that they would still 
hold, mutas mutandis, after replacing everywhere the retarded prescription with an advanced 
one: The corresponding classical field A % would have support only at negative x~ , while the 
semi-fast gluon exchanged within x would freely propagate at positive x~ . 

In Refs. [||,[14],[il| , where the advanced prescription has been used extensively, one has 
found similar simplifications in the calculation of Feynman graphs for, e.g., the scattering 
of a quark or a gluon off a hadronic target. As shown there, it is only with this prescription 
that one can ignore the final state interactions of the struck quark (or gluon) in deep inelastic 
scattering, which is essential if the produced jet is to be used as an indicator of the hadron 
wavefunction. Our present analysis corroborates the conclusions in Refs. 0,[14|,[16| that 
retarded and advanced LC-gauge prescriptions are special in that they not only lead to 
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technical simplifications, but also allow for a more transparent physical interpretation of 
the results. With more symmetrical prescriptions like principal value PV^- or Leibbrandt- 
Mandelstam, one cannot avoid the overlap between the classical fields and the quantum 
fluctuations, and thus neither the final/initial state interactions. 

To summarize, when computing \ with a retarded prescription, one can replace G^ v 
in eq. Q2.1|) by the corresponding free propagator Gq V . This calculation will be presented in 
the following subsection. For completness, in Appendix B we shall verify, on the example of 
X3, that a lengthier calculation using the full propagator leads eventually to the same 
result. 



B. Explicit calculation of x{x,y) 

With G l i replaced by the free propagator Gq = S^Gq the first term xi i n the r.h.s. of 
eq. (|2.1|) has been already computed in Sect. 5.2 of Paper I, with the following result: 

Xi = £ In(l/6) Ft j J |^ . (2.7) 



Consider now the matrix element (|2.2j ) which enters X3- With G ++ — > Gq + , cf. eq. ( |A10|) , 
this gives (for y + = x + , and x~ ~ y~ ~ 0) 

io + + te io + — te 

strip 27T J 27T J {2^ P + + 16 JT ^ P + ~ 16 ' ^ 

where the integral over p~ is restricted to the strip (|1.31| ). The various factors of p + in the 
integrand can be combined into a PV prescription in l/p + : 

2p+ 11 1 . . 

I — = 2PV — . (2.9) 



(p + + it) (p + — ie) p + + ie p + — ie p 
Then, the integral over p + is easily computed by contour techniques: 

* + FvJL L, = ^>, (2. 10 ) 



27r p + 2p + p~ — p 2 ± +ie 2p\ 

where e(p~) = 0(p~) — 0(—p~) is the sign function. The restricted integral over p~ generates 
the expected logarithmic enhancement (below, A = A - , cf. eq. ( |1.31|) ) : 

dp' e(p~) ( r- A r A / b \ dp' e(p~) 1 , , ,,, , 

+ / ^— = - ln(V & )- ( 2 - n ) 



/strip 27T p~ JA ) 2n p~ 7T 

The final result reads 
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(0,x± 



-J~-G ++ -^—%y±) = -ln(l/6) / ^%\^ x ^ 
io + + xe io + — le it J (2-k) 2 p z ± 



(2.12) 



which immediately implies: 

o 2 f d 2 p^_ 1 

Xs = — ln(l/6) / 777-7^ — 

A similar calculation yields: 



(2.13) 



<0, acjL| 



G + %y ± ) = (0,x ± \G 



id + + ie 



id + — ie 



\o,y±) 



(2.14) 



which allows us to also compute X2 
- 2 - d 2 P± 1 



X2 



7T 



ln(l/6) 



(2tt)2 pi 



(2^2^ + p) x e 



»pj_-(sj_-y_L) 



(2^>+^') 3/ + 



+ {2F +i d l ) x jP^-v^ (2V jj F +j + p) y \ (2.15) 



The previous results for Xh X2, and X3 are conveniently combined as: 



X(S, $ = ^ hi ~ [r? 6f(x ± - y ± ) T+i + 



71 b 



2T +l \V i )+p 



(x±\—- T \y±) 



2(V^-^)^ + p 



, (2.16) 



which is our final expression for x- The following notations have been used: 

d 2 p± 



$±(x± - y± 
1 



[x± 



-V 2 



y±) 



(2vr) 2 v 
d 2 p ± 1 
(2?Z 



5 lJ 



p l p> 

pI _ 



3 ipx'(a:±-2/_L) 



jp±-(x±-y±) 



(2.17) 
(2.18) 



C. The two-dimensional correlator x( x ±,y±) 



In eq. ( [2.161 ), both p and have support near the LC, at ^ x~ = 1/A + . 
Thus, although induced by quantum modes with relatively large longitudinal wavelengths 
Ax~ ^> 1/A + , the charge-charge correlator appears to be as localized in the longitudinal 
direction as the original source at the scale A + . This is so because the vertices responsible 
for this quantum effect are explicitly proportional to p or JF +J (cf. eq. (|1.26| )). Thus, x^x.y) 
is manifestly sensitive to the internal structure of the source, and also to the structure of the 
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background field at small distances x ^ x T . (Note, in particular, that eq. ( |2.16|) involves 
the product J r+l (x)A l (x), and thus the field A l (x) within the support of p.) 

If this was also true for the two-dimensional density x( x ±> V±)i which is the quantity 
which enters the RGE ( |1.16|) , this would spoil the separation of scales assumed by the 
effective theory, and thus the validity of the latter. Note that, the simple fact that x{ x ±-> V±) 
is obtained from after integrating out x~ and y~, cf. eq. ( |1.32| ), is by itself not 

sufficient to guarantee that the physical information about small x~ or y~ is truly irrelevant. 
For instance, the following integrated quantity 

J dx-F +l (x)A\x), (2.19) 

although a function of x± alone, is nevertheless sensitive to the values of A l (x) at x~ ^ x~, 
since the electric field T +l (x) has its support there (cf. eq. ( |1.7| )). On the other hand, the 
following quantity (recall that A l {x~ — > oo) = Alo, while A l (x~ — > — oo) = 0) 

J dx~F +i (x) = J dx~d + A l = A i 00 (x ± ) (2.20) 

is not sensitive to the structure of A l (x) around the origin, but only to its asymptotic value 
at large x~ . 

As we shall demonstrate now, this latter situation applies also to x(x, y), which, like 
the electric field total derivative with respect to its longitudinal arguments x 

and y~: 

X (x,y)=d+d+ C(x,y) . (2.21) 

To see this, note the following identities, which hold for an arbitrary function $(xjj, 

{2F H V i + p) x $(x ± ) = id+Vmx L ), 

S(yj_)(2P^ + p) y = -t<S>(y ± )Vi 2 d+ (2.22) 

(with dy in the second line acting on the function on its left). These identities rely on the 
classical equation of motion (\1.3j) , that is, 

&F +i -ig[A i ,J rH ]=p. (2.23) 

For instance, one can obtain the first identity by writing: 

d + V 2 ^{ Xi _) = 9+ ((9* - igA l ){d% - igA'ty) 

= - - ig{T H A i + A* J**)® 

= -i{d l F +l + %g[T Vl A l - A l T+ l ) + 2^ +l D i }$ 
= -i(p + 2F +i V l )<$> (2.24) 

where in writing the second line we have used <9 + $ = and d + A l = J-" +t , and the last line 
follows from the previous one after using Q2.23 ). 
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By using eqs. (|2.22|) and ( |2.16|) , one deduces that x{x,y) is indeed of the total- 
derivative form ( [2.21|) , with 

C{x,y) = ^lni 6f(x ± -y ± ) A{ + (P'U 4 ),^!^!^)^^'),} • (2-25) 



Then, the integrations over x and y in eq. (|1.32| ) become trivial, and yield 



x(x±,v±) = - y±)^L(y±) + (vlAl) x Ax±\^\y±)(Alv%) y± } (2.26) 

(with = d % — igA* oc (x±)) which is sensitive only to the asymptotic fields. 

The gauge rotation ( |1.21|) of eq. ( |2.26|) to the background COV-gauge is straightfor- 
ward. By using A z oa (x±) = [if ' g)V x d l V^ (cf. eq. Q1-6Q ), one immediately finds 

<r\:^{.r. -y ± )<PV y + di((dW)jx ± \^\y ± )(&V) y yj}, (2.27) 



x{x±,y±) = — 
9 

where = V^(x±), cf. eq. (|1.8Q , and the derivatives not included in the brackets act on all 
the functions on their right (or left). 



D. From \ to rj 

Eq. fl2.27| ) provides x(xj_, y±) as an explicit functional of the field a(x), via the asymp- 
totic Wilson lines V(x±) and V*(x±). It is therefore more convenient to write down the 
evolution equation in terms of a rather than p. The appropriate transformation x ~~ * V is 
shown in eq. (|1.23b ), and will be worked out in detail in what follows. 



We shall do that in two steps, corresponding to the two terms in the r.h.s. of eq. ( [2. 27 ) 



which we write as x — Xi + Xi- Correspondingly, r] = r\\ + r^ 2 with (we omit the factor of 
4/g 2 at intermediate steps) 

m(xx,V±) = Jd 2 z L JdW(x ± \^\z ± }(dldi{v}6Y(z ± -u ± )V u })(u ± \-^\y ± } 

d 2 z ± J d 2 u ± (x ± \-^\z ± ) Vl (z ± \5 l > - -^g- \u ± ) V u ( u± \-^\y ± ), (2.28) 

where the second line follows after some integrations by parts, and we have written 



Similarly, 
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ife(*L,V±) = - [d 2 z ± fdW(x ± \^\z ± }(d^)(z ± \-^\u ± )(d^V u )(u ± \^\y ± } 

d l 1 , d j 



By using 



d 2 z x J d^V^dldiix^-^lz^iz^-^MM-^ly^. (2.30) 



di(x ± \ \y ± ) = -8®(x ± - y ± ), (2.31) 



this can be further transformed as 
V2(x±,y±) = Vl (x ± \ —^2- |yj_> V y 



+ J d 2 u ± VI <rrj_| — \u±) V u (u ± \ — ^- |y ± ) 

+ / d% ± J d 2 z x (x x \ \z x ) Vl (z ± \ ^ \u ± ) KM |yj • (2.32) 



The last term in the above expression is non-local in both u and z, but it cancels against 
the similar term that contributes to rjx, eq. (|2.28 ). So we arrive at 



v(x±,y±) = (i + v^v y ) (x ± \ l^> 



+ / d 8 *!. <arj.| (z ± \ ^ ly±> (^K + V?y») • (2.33) 



By also using 



<l 2 z± (x ± \ [«jl> (z±\ \y±) = -(x ± \ \y ± ) (2.34) 



we finally express r] as (below, we reintroduce the overall factor 4/g 2 ) 



The kernel in this equation has been written in coordinate space by using eq. ( p.29| ) and 



d 2 P± f d 2 k ± P± -k ± e i p± .( Xx - Zx ) e i kx .( Zx -y x) = 1 & - Z*)^ - Z*) (2 3g) 



(2tt) 2 J (2vr) 2 p\k\ (2tt) 2 (x ± - z ± ) 2 (y ± - z ± ) 2 

Eq. (|2.35 ) is our final result for the "real correction" , to which we shall return in Sect. |IV A 
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III. THE INDUCED COLOR CHARGE, OF "VIRTUAL CORRECTION 



We shall now similarly compute the coefficient v in the RGE (|1.22|), or the "virtual 



correction". This calculation will be organized as follows: By evaluating the r.h.s. of 
eq. fll.36 ), we shall first obtain the three-dimensional color charge density &(x) in the LC 



gauge. The longitudinal structure of the result will be explored, and then integrated out, 
to yield the two-dimensional density a a (x±), cf. eq. (|1.35|) . In a separate calculation, based 



on eq. fl2.27|) for we shall obtain the classical polarization Sa(x±) defined by 

eq. ( |1.21c|) . By combining the various results as shown in eq. ( |1.21b|) , we shall deduce the 
induced source cr a (x±) in the (background-field) COV-gauge, from which v a (x±) will be 
finally extracted, as shown in eq. (|1.23a| ). 



At intermediate steps in this calculation we shall meet with logarithmic "tadpoles", 
i.e., contributions to a or 5a which are local in x± and proportional to 

(3-D 

P± 

Moreover, the momentum integrals giving the remaining, non-local, contributions will gener- 
ally involve logarithmic ultraviolet (UV) divergences. After summing up the various contri- 
butions, all UV divergences must go away. Indeed, our one-loop calculation cannot introduce 
other UV divergences than those of the standard perturbation theory, which, moreover, must 
be absent as well, since they are not enhanced by the longitudinal logarithm \n(l/x). (In 
other terms, the UV divergent terms are purely 0(a s ) effects, which are a priori ignored in 
our LLA.) To shorten the presentation, in this section we shall simply ignore all the tadpoles. 
At the very end, we shall verify that the final, non-local, result for a a (or u a ) is UV finite, 
meaning that all the tadpoles must have cancelled among themselves. 

The infrared finiteness, on the other hand, is a different issue, which, as in the cor- 
responding BFKL problem, is related to subtle compensations between real and virtual 
corrections. We shall return to this issue in the discussion of the results in Sects. [TV] and M. 



A. Explicit calculation of a(x) 



Eq. ( |1.36| ) shows that there are two contributions to the induced source: a = &i + a 2 . 
The second contribution <j 2 , to be computed in Appendix C, turns out to be a pure tadpole, 
and as such it will be ignored in what follows. The first contribution &i, to be computed 
below, is given by: 



Tr(T a ax(f)), a^x) = -gd+G l \x } y) 



(3.2) 



It can be verified that, out of the four different terms in eq. (|A12a|) for G u , it is only the 
last one, namely 
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G»(x,y) = -(x\V^G ++ ^V^\y) (3.3) 

/■ 1 1 

= Z>* / rizfrfz^ (aT | [zr)G ++ (zr, z ± ; z 2 ~, y ± ) (z 2 ~ | _ |y~)X>^ , 

that gives a non- vanishing contribution to &\. The derivative dy in eq. (|3.2|) is readily 
performed as : 

#terl^:lv-> (sj* + ^($) = -y-)2?t* + ^ 2 1^l<r>^ + X£), (3.4) 

thus giving the following two contributions to <j\ : 

a x (x) = -ig^x\V' l -^G ++ -^\x)F + \x) - ig{x\V i j^G ++ V*\x) = a n + a 12 . (3.5) 



The first term <7n involves the same matrix element ( |2.4j ) as %3 i n ec L- ( P ■ 1| ) • Indeed, since 
jF+'(a;) is localized at x~ ^ x~, one can effectively set x~ ~ in the above expression <7n, 
which leads us to the matrix element evaluated in eq. ( 12.121) . This gives 



n 2 r d 2 r> < piPx-fax-vx) 
hi{x) = - 9 -Hl/b)Vl / ^ e —- 2 



7T v ' ' x J (2tt) 2 p\ 

d 2 p ± 1 



x=y 



= i 9 - Hl/b)A\x)^\x) Jj^y 2 ^, (3-6) 

which turns out to be a pure tadpole, and will be therefore ignored in what follows. As for 
the second term <t 12 in eq. (1331), that is, 
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-ig(x\V— -—G ++ V^\x), (3.7) 
io + + it 



this is determined exclusively by the "crossing" piece G ++ ^ of the propagator G ++ , as 
given by eq. (|A8|) . (The corresponding "non-crossing" contribution is shown to vanish in 
Appendix C.) 

To summarize, the whole non-trivial contribution to a a {x) comes from the crossing 
piece of eq. (|3.7|), which after using eq. ( |A8| ) can be rewritten as (with d\ 2 simply denoted 
as a from now on) 



a = -ig 



dp- 2i r dp+ dk + e -^ p+ - k+ > 



strip 2tt p- J 2tt 2tt p + + ie 
xVlJ dT ± ( P± -k ± )G (p)G (k) 

x {6(p-)V(x ± )V\z ± ) -6(-p-)V(z ± )V\y ± )\vl l , (3.8) 



where it is understood that k = p , and we have used the following shorthand notation 
for the transverse phase-space integrals (this is a function of x± and y±): 
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dT, = /fe /^LL /dV e 



(2vr) 2 7 (2tt) 



(3.9) 



Thus, the transverse derivatives V l x and in eq. fl3.8|) act on both G?r_i_ and the Wilson 
lines, and the limit — >■ x± is taken at the very end. 



The integrals over p + and k + in eq. (|3.8| ) are readily performed via contour techniques: 



dp + 
~2V 



1 



P± 



p + + ie 2p + p~ — pj_ + ie 
x~) - \e(x~)6(p-) - 6(-x-)6(-p-) 



e 2 p~ 



(3.10) 



where the first term 9{x ) in the accolades is the contribution of the LC-gauge pole at 
p + = 0, while the other terms comes from the single-particle pole at 2p + p~ = p\. Similarly: 



dk+ 
~2^ 



1 



2k + p- -k± + ie 2p~ 



[e(p~)6{-x~) - 6{-p-)9(x~)) 



(3.11) 



In the product of eqs. ( |3.10|) and ( |3.11|) , the only term which survives is the one proportional 
to 9(—p~)9(x~) survives. It yields 



&(x) = g6{x-)Vi 



dp 



-p~ 



strip 2ll (p 



-\2 



dT 



p ± ■ k± i±L_ 



p 2 ± 



e**- x ~V{z ± )V^{y ± )Vl i 



y=x 



(3.12) 



with the important consequence that the induced source has support only at positive x~ . 
Clearly, this property is related to our retarded gauge-fixing prescription, which is responsible 
for the 9{x~) term in the r.h.s. of eq. ( |3.10| ). If an advanced prescription l/(p + — ie) was 
used instead, eq. (|3.10|) would have rather generated a term in 9(— x~), and thus a quantum 
correction &{x) with support at negative x~ . 

To further explore the longitudinal structure of the induced source, it is convenient to 
perform the restricted integration over p~ in eq. Q3.12 ). One thus obtains: 



dp 



strip 2ir (p 



p ) 

e 2 p 



1 e ~ ibA+x -e 



iA+x~ 



-\2 



mx 



(3.13) 



where we have replaced k\/2K « A + to LLA (cf. the discussion after eq. ( 1.31 )). This 
shows that the ^"-dependence of a is carried by the following "form factor" 



-ibA+2 



F(x' 



9(x~) 



— e 



-iA+3 



X 



(3.14) 



which has support at 1/A + ^ x~ ^ 1/6A + . Indeed, F(x~) ~ both for small x~ <C 1/A + 
(since in this case the two exponentials mutually cancel), and for large x~ ^> 1/6A + (where 
the two exponentials are individually small). For x~ ^ 1/A + , T> % = P^, with (cf. eq. (|1.6| )) 
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= & - igdLM = 9* + Vd { Vl 



(3.15) 



By also using 



(e- ik ^ y± V^(y ± ))v^ y = (d i e- ik± - y± )V" ! (y ± ) = ik i e^ ik± - v± V i (y ± ), 
eq. ( |3 . 1 2| ) finally becomes (with T> % = V^) 



a(x) = 9 -F{x~)Vi J dTj-^VV z V} 



(3.16) 



(3.17) 



which shows explicitly that the induced source cr{x) has support at 1/A + ^ x ^ 1/6A + , 
as anticipated in eq. ( |1.18| ). Moreover, the logarithmic enhancement is not yet manifest on 



eq. ( 3.17 ). This comes out only after the integration over x , via the following equation 



dx Fix 



'4 



(3.18) 



which together with eq. (|1.35| ) implies: 

ga(x ± ) = 4 ( V 



dv , p± ' k± p v Wt 



(3.19) 



In this equation, the limit y± — > x± within dT± can be taken also before acting with the 
derivative V> % x . Indeed, it can be easily checked that the additional term generated in this 
way cancels out after taking the color trace with T a . Thus, in eq. ( |3.19|) , 



dT± = J (2^ J (2^ J d2z±e 



(3.20) 



For what follows, it is useful to have the expression of a rotated to the covariant gauge. By 
using V^V l O = d % (V^O), one obtains: 



g(V^V)(x ± ) 



491 1 dT± l^t kiV * Vz - 



(3.21) 



Eqs. ( |3.19| ) and ( |3.21| ) represent our final results in this subsection. But for what follows, 
it is important to also keep in mind the longitudinal structure displayed in eqs. ( |3.14 ) and 
( p,17| ). These equations show that, unlike the real correction x{x,y) which is as localized in 
the longitudinal direction as the original source p(x), the virtual correction a(x) is relatively 
delocalized, and has support on top of the original source. As discussed in Sect. [TB| , this 
justifies the peculiar longitudinal structure of the RGE (|1.16|) or (|1.22|) . 
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B. From a to v 



According to eq. ( |1.21b[ ), the induced source a a (x±) in the (background- field) COV- 
jauge is obtained by combining 



Vi(x ± )a b (x ± ) = Tr(T a V} a(x ± )V x 



(3.22) 



with the classical polarization term defined by eq. ( |1.21cj ), that is, 

5a a {x ± ) = Tr{T a 5a{x ± )) 

Of 1 

Sa(x±) = ?- j d y± x(x±,y±) (y±\ — ^- \x±) = 5ax + Sa 2 



(3.23) 



with the two pieces 5a\ and 5a 2 corresponding to the two terms in the r.h.s. of eq. (|2.27|) 
for x- Their calculation is very similar to that of rji and i] 2 in Sect. [11 0. The results are 



U<^ l U ± ) = 2 j dl^^^p^V^V, 



g5a 2 (x±) = 2 dT 



pj_k' 
(pi - k ± ) ■ k ± 



p 2 ± k 2 



{p\d l V%V z + i(V 2 ± V%V z } , (3.24) 



with dT± defined as in eq. (|3.2(J|). After adding these two expressions together, and using 



/ dT ± (?A#vi) x v z + i(v 2 ± v%v 2 ) = id x J dr ± V A^- (d^) x v z , (3.25) 



P± ■ kj 



pj_k 



p{kj 



(which follows from the symmetry of the integrand under the exchange p l 
obtains: 

gSa(x ± ) = 2\i<f x [ d r ± ^±(d i V%V z - i(V 2 ± V%V x 1 <h>L ' 



p\k] 



(27T) 2 pi f ' 



-k l ) one 
(3.26) 



where the second term within the braces is a pure tadpole and thus will be neglected here, 
to keep in line with our previous strategy. Our final result for the classical polarization term 
is therefore: 



g8<T a {xx) = 2idl I dT ± ^ r ^d i x Ti(T a V^V z 



(3.27) 



As a matter of fact, the transverse integrations in eq. ( |3.27|) develop a logarithmic ultraviolet 
divergence, which in the original expression (|3.26|) was precisely cancelled by the explicit 
tadpole there. Still, neglecting this tadpole is a legitimate procedure since the UV divergence 
in eq. ( |3.27| ) is actually needed to cancel a corresponding divergence in eq. ( |3.21| ) for V^aV 
(see below). This simply shows that in a more complete calculation where all the tadpoles 
would have been kept explicitly, the tadpole in eq. ( |3.26| ) would have cancelled against 
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similar tadpoles emerging in the calculation of a, and which have been neglected in arriving 
at eq. flOlp . 



To conclude, the final result for the induced source a a (x±) in the COV-gauge is ob- 
tained by combining eqs. ( |3.21|) and ( |3.27| ) according to eq. ( |1.21b|) , and reads: 



gt a (x ± ) = -V 2 ± [2i J dT ± ?^Tr(T a VX)} = -gV 2 ± v a (x ± ) 



(3.28) 



In order to obtain this simple form, it has been convenient to rearrange the integrand in 
eq. (|3.21|) as follows: 



dT ± ---k i = (1/2) dT ± --- (k 1 - p 1 ) = H/2R / dT ± 



(3.29) 



which relies again on the symmetry of the integrand under p % <-> —k l . 

As anticipated, the expression ( |3.28| ) is UV finite. This can be verified as follows: Let 
q± = p± — kj_ be the momentum conjugate to the transverse distance x± — z±. Then 



P±kj 



X 



2bi 



(pj 



pKpa 



where the following identity has been used: 



P± ■ {Pj 



pKpa 



q1 



(Pa 



pKpa 



(3.30) 



(3.31) 



Since q± is kept finite by the non- locality in x± — z±, potential UV problems can occur 
only from the limit p± — > oo at fixed q±. In this limit, the last term within the braces in 
eq. (|3.30|) is manifestly UV finite (since it decreases like l/p±), while the first two terms 
give logarithmic tadpoles / d 2 p±/p 2 L . But in these two terms, the integration over q± sets 
x± = z±, so that VjV z — > 1. Thus, the potential tadpoles go away after taking the color 



trace with T a . Incidentally, the above discussion shows that eq. fl3.28|) can be equivalently 
rewritten as: 

d 2 p± f d 2 q ± 



9<r a {x± 



d'z ± 



Jq±-(x±-z±) 



(2vr) 2 "* Pl(p±-q±) 2 " 



(3.32) 



(2vr) 2 

a form which will be useful later on. 

From eq. ( 13.28 ), the coefficient z/ a (x_i_) in the COV-gauge RGE ( |1.22| ) can be trivially 
identified, according to eq. (|1.23a|) . It reads: 

d 2 z± 



gv a {x ± ) = 2i 



(2tt) 2 (x ± -z ± )' 



Tr(T a VjV z 



(3.33) 



where eq. Q2.36| ) has also been used. Physically, the quantity is a (x±) is the induced color 
field in the COV-gauge, that is, the modification in the original field a a induced by quantum 
corrections |L1 . 
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IV. THE RENORMALIZATION GROUP EQUATION 



We are finally in a position to write down the renormalization group equation explicitly 
and study some of its properties. In the a-representation, this is given by eq. ( |1.22| ) with 
the coefficients r\ and v from eqs. ( |2.35| ) and ( |3.33| ), respectively. We summarize here these 
equations for convenience: 



dW r [a] 
dr 



2 5a^(x±)Sa b T {y A 



-[W T 7] ab 



xyl 



5a°:(x± 



■M 



(4.1) 



where 



d 2 z± 



x 



z i )(y i 



2tt J 



7r J (2tt) 2 (x± - z±) 2 (y± 
d 2 z ± 1 



.{i + vX-v ] x v z -v}v y } a \ 



(2tt) 2 ( x± - z± y 



Ti[T a V^V z ). 



(4.2) 
(4.3) 



Note that, as compared with the original equations (|1.22| ), ( |2.35| ) and ( |3.33|) , we have 
rescalled here the coefficients r\ and v to absorb the overall factor a s = g 1 j^is. The equations 
above represent our main result in this and the accompanying paper |[L1|| . They govern the 
flow with t = ln(l/a;) of the probability density W T [a] for the stochastic color field a a {x) 
which describes the CGC in the COV-gauge. We now turn to a systematic discussion of the 
properties of these equations. 



A. General properties 



i) The coefficients rj and v are real quantities. This is so since the Wilson lines in the 
adjoint representation are real color matrices: V* b = V a b, and therefore V^ h = Vb a - By using 
this same property, one can further verify that: 

ii) The 2-point function r] ab (x±,y±) is symmetric and positive semi-definite: 

>f(ZL,</±)=>J fc, (</±,I±) 

7T J (2tt) 2 (x± - z±y{y± - z±) 2 v > v ' 

This guarantees that the solution W r [a] is positive semi-definite (as it should) at any r 
provided it was like that at the initial "time" To- 

iii) The RGE preserves the correct normalization of the weight function: 

J Va W T [a] = 1 at any r. (4.5) 



Indeed, the r.h.s. of eq. ( |4.1| ) is a total derivative with respect to a. Thus, if eq. ( |4.5| ) is 
satisfied by the initial condition at time Tq, then it will be automatically true at any r > r . 
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iv) The solution W T [a] encodes the information about the longitudinal support of 

a a (x~,x±) and its evolution with r. From the calculation of quantum corrections, we know 
that the quantum evolution up to r generates a non-trivial color field a a (x) only within 
the range < x~ < x~ , with x~ = XqQ t and Xq = 1/P + . Thus, as a functional of 
{a a (x~ ,x±) | < x~ < oo}, W T [a] must have the following structure: 

W T [a] = 5[ a> ] W T [a<]. (4.6) 

Here, a< (a>) is the function a(x) for x~ < x~ (x~ > x~), 

a(x) = 9(x~ — x~)a < (x) + 9(x~ — a;~)a>(x), (4.7) 

the functional 5-function 5[a; > ] should be understood in terms of some discretization of the 
longitudinal axis, as the product of ordinary (^-functions at all the points x~ > x~ : 

5[ a> ] = I] <*(<**-) (4-8) 

(this simply shows that a a (x~ ,x±) = with probability one for any x~ > x~), and the new 
functional W r [a<] involves a a (x~ ,x±) only within the restricted range < x~ < x~. 

The structure ( |4.6|) can be used to simplify the average over a when computing ob- 
servables: 

(0[a]) T = f VaO[a]W T [a] 

J0<x-<oo 

= f VaO[a]W T [a], (4.9) 

J0<x <x T 



where in the first line the functional integral runs over fields with support at x~ > 0, while 
in the second line the support is restricted to x~ < x~ . The expression in the first line is 
nevertheless more convenient to derive an evolution equation for (0[a\ ) T , since there the 
whole dependence on r is carried by the weight function: 

§p(0[a]) T = J VaO[a]^^, (4.10) 

with dWr/dr determined by the RGE Q). 

v) Within the RGE, the Wilson lines and V, which were a priori defined as path-ordered 
integrals along the whole x~ axis (cf. eq. (|TT%|)), can be effectively replaced by: 

V\x L ) -> U j T {x ± ), Ul(x±) = Pexpj^^ dx~ a{x~,x ± )^j . (4.11) 

This obvious consequence of eq. fl4.6p is interesting since it shows that the functional deriva- 
tives within the RGE ( fll|) act on Wilson lines as derivatives with respect to the color field 
at the end point x~ —x~. Explicitly: 
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gg = ^ - ,j ™t (l±) , = _ 19 „, (s± _ z±) UMT . (4 . 12) 

which should be compared with more general formulae like: 

= ^(* ± - z ± ) Ul !T (x ± )T a Ut(x ± ) , (4.13) 

where the integral over x~ within T runs from x~ up to oo. Clearly, eq. ([4.13 ) reduces 



to the first eq. ( 4.12 ) for a field a with support at a; < x T . 

vi) The longitudinal coordinate and the evolution time get identified by the quantum 
evolution. The previous arguments at points iv) - vii) reveal a strong correlation between 
the longitudinal coordinate x~ and the rapidity r which eventually allows us to identify 
these two variables. To make this identification more precise, let us introduce, in addition 
to the momentum rapidity r = ln(P + /A + ) = ln(l/x), also the space-time rapidity y , defined 
as (for positive x~ and Xq = 1/P + ): 



Then the field 



y = ln(x /x Q ) , — oo < y < oo . (4-14) 



®y(x±) = a a (x =x y ,x ± ), x y = x e y , (4.15) 



at positive^ y (i.e., at x~ > Xq ) has been generated by the quantum evolution from r' = 
up to r, with a one-to-one correspondence between y and r. That is, the field ay within 
the space-time rapidity layer y < y' < y + Ay has been obtained by integrating out the 
quantum modes within the momentum rapidity layer r < r' < r + Ar with r = y and 
At = Ay. This leads us to treat y and r on the same footing, as an "evolution time". 
Then {a y (x±) \ — oo < y < oo} is conveniently interpreted as a trajectory in the functional 
space spanned by the two-dimensional color fields a a (x±) ; this trajectory effectively ends 
up at y = r. With this interpretation, eq. ( |4.1|) describes effectively a 2+1 dimensional field 
theory (the two transverse coordinates plus the evolution time), which is however non-local 
in both x± and y (since Wilson lines likes fl4.11|) involve integrals over all y G (— oo, r)). 



vii) The infrared and ultraviolet behaviours of the RGE. These are determined by the 
corresponding behaviours of the kernel t] ab (x±,y±) (we shall shortly see that v is related to 
77), and, more precisely, by the behaviour of the integrand in eq. Q4.2| ) for fixed x± and y± 
(since these are the "external" points at which we probe correlations in the system). For 
the infrared, we need this integrand in the limit where z± is much larger than both x± and 
y±. Then the products of Wilson lines involving z± are expected to be small (since, e.g., 
(V^V Z ) — > as \z± — x±\ — ► 00 f2~l"H), so it is enough to study the large-zj_ behaviour of 



8 The field a y at negative y, i.e., at x < x , is rather a part of the initial conditions, in that it 
exists independently of the quantum evolution. 
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lC(x ± ,y ± ,z ± ) = Z [lf Z%) 2 . (4.16) 

For z± ^> x±, y±, this gives JCir ~ 1/^1 and the ensuing integral (<i 2 ;z_i_/r[) has a logarithmic 
infrared divergence. Thus, there is potentially an IR problem in the RGE. This is not 
necessarily a serious difficulty, since IR problems are expected to be absent only for the 
gauge-invariant observables. In this context, we would expect the IR divergences to cancel 
out when the RGE is used to derive evolution equations for gauge-invariant observables. 
Although a general proof in this sense is still missing, we shall nevertheless see some explicit 
examples where such cancellations take place indeed. 



Coming now to the ultraviolet, or short-range, behaviour, it is easy to see on eqs. ( |4.2|) 
or Q4.4Q that no UV problem is to be anticipated. For instance, the would-be linear pole of 
fc(x±, y±, z±) at \z± — x±\ — > is actually cancelled by the factor 1 — V£V X which vanishes 
in the same limit. 



B. The RGE in Hamiltonian form 



The RGE fl4.1|) is a functional Fokker-Planck equation, that is, a diffusion equation for 
a (functional) probability density. It portrays the quantum evolution towards small x as a 
random walk (in time r) in the functional space spanned by the two-dimensional color fields 
a a (x±). The term involving second order derivatives in this equation is generally interpreted 
as the "diffusion term" , while that with a single derivative is rather the "force term" . But 
this interpretation is not unique for our RGE where the 2-point function r] ab (x±, y±) — the 
analog of the diffusion "constant" — is itself a functional of the field variable a a (x±), so 
that eq. (|4.1|) may be as well rewritten as 



= — i — |i «» _sw L _ i\ — \ i 

which rather features ^(Srj/Sa T ) — v as the effective "force term". An important property, 
which is full of consequences, is that this effective "force term" is precisely zero. 

Specifically, the following functional relation holds between the coefficients of the RGE: 

It is straightforward to prove this relation by computing the effect of 5/5a b T {y±) on the terms 
involving Wilson lines in eq. ( J4.2[ ) for r) ab (x±,y±). For instance, 

r) fiV^ ac SV cb 



Sa b (y ± ) " x v > 6a b (y ± ) » - Softy ± ) 

= ^g^\x 1 _ - y x )(T b V})j? - ig5^(0 ± )V^ c (y y T b ) = , (4.19) 
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where each of the two terms vanishes because of the antisymmetry of the color group gen- 
erators in the adjoint representation (e.g., (T b ) a {, = 0). One can similarly show that: 



8 



Sa b (y ± ) 
5 



(Vx V z 



ab 



-igt 



5 (2) (x ± -y ± )(T b V2V z ) ab = ig5^(x ± -y ± )Tr(T a V^V i 



(y}v y ) ab = o. 



(4.20) 



The only non- vanishing contribution is that in the first line above, and this precisely repro- 
duces eq. (|4.3|) after integration over y±, since (cf. eq. (|4.16|) ) 



JC(x ± ,x±,z_ L ) 



1 



(4.21) 



(x± - z ± ) 2 

By using eqs. ( 4.17 ) and (|4.18 ), the RGE can be finally brought into the Hamiltonian form: 

dW T [a] 



Or 



HWJa] 



(4.22) 



with the following Hamiltonian: 

1 S S 

JiM - / ^ Y~~~~~~To (i - vxu 



d 2 z ± ■ 



2tt 



r a (z±)r a (z ± ), 



2/T (z± — x±) 



Sa b T (x ± ) 



(4.23) 



which is Hermitian (since r^ y is real and symmetric) and positive semi-definite (since the 
"current" J l a {z±) is itself Hermitian). 

At a first sight, the interpretation of eqs. ([4.22 )-( 4.23 ) as a Hamiltonian system may 
look compromised by the fact that the operator ( [4.23| ) appears to be time-dependent, e.g., 
via the T-dependence of the differentiation point a®(x±), and even non-local in time, via 
the Wilson lines V and V* (which involve integrals over all times y < r; cf. the discussion 
after eq. ( |4.15|) ). But this rather means that we have not yet identified the correct canonical 
variables: these are not the fields a", but rather the Wilson lines ending at r. Specifically, 
the canonically conjugate variables and momenta are {U® b (x±), U^(y±)}, with U T (x±) as 
defined in eq. (|4.11| ), and n"(xj_) acting on U T or U\ as the functional derivative w.r.t. the 
field at the end point: 



KM = - * 



ig 5a"(x±) 



(4.24) 



The associated Poisson brackets, or equal-time commutation relations, are easily inferred 
from eq. (|4.12 ): 



n«(x ± ), U}(y ± )] = T a Ul(y ± )5 (2 \x ± - y ± ), \u a T (x ± ), U T (y ± )} = -U T (y ± )T a 5^ (x ± - y ± ), 



U a T (x ± ), U b T (y ± )] = if ohc U c T (x ± )6^(x ± - y ± ), 



(4.25) 
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(Non-shown commutators like [U$, U±] or [UJ, U T ] vanish as well. This should not be confused 
with the fact that, as color matrices, U T (x±) and U T (y±) do not commute with each other 
at different points x± ^ y±.) 

Eqs. (|4.11 ) and ([4.24]) provide explicit representations for the canonical variables in 



terms of the gauge field a^(x±). But we can extend these relations to more abstract defini- 
tions, which hold independently of any representation. The generic canonical "coordinates" 
are the group- valued 2-dimensional fields V (x±), while the canonical "momenta" n a (xj_) 
are Lie derivatives w.r.t. these fields, whose action is defined by eqs. (|4.25|) . These momenta 
satisfy the commutation relation of the Lie algebra, as they should (cf. the second line of 
eq. ( |4.25| )). In terms of these variables, the Hamiltonian ( |4.23| ) takes the more standard 
form of a second-quantized Hamiltonian for a 2-dimensional field theory: 

H[U, V] = \J d 2 x ± J d 2 y ± U a (x ± ) (g 2 V ? y [V, 0]) n V) • (4-26) 

Note that this is a purely kinetic Hamiltonian (no potential). It describes free motion on a 
non-trivial manifold (here, a functional group manifold), with the kernel fi^^fV, V>] playing 
the role of the metric on that manifold. Then, (|4.22j) is like the evolution equation for this 
Hamiltonian in imaginary time. Actually, giving the probabilistic interpretation of W T [a], 
the best analogy is not with the (functional) Schrodinger equation in imaginary time, but 
rather with the Fokker-Planck equation. We shall return to this analogy by the end of this 
section (see also Ell). 



The Hamiltonian ( |4.26| ) has been first obtained by Weigert [20| in a rather different 
context, namely from an analysis of Balitsky's equations. These are coupled evolution 
equations for Wilson-line operators that have been derived by Balitsky |TS| via the operator 
product expansion of high-energy QCD scattering in the target rest frame. Weigert has 
subsequently recognized that all these equations, which form an infinite hierarchy, can be 



generated from a rather simple functional evolution equation with the Hamiltonian ( (4.26| ). 



The fact that we have come across the same Hamiltonian means that our RGE too generates 
Balitsky's equations, which we shall verify explicitly in the next section on the example of 
the 2-point function (tr(V^ (x±)V(y±))) T . 



Weigert has been also the first one to notice the remarkable relation ( J4.18 ) between 



the coefficients in the functional evolution equation. Let us devote the end of this section 
to a more thourough discussion of this relation and some of its consequences. 



As a relation between one- and two-point functions, eq. ( 4.18 ) is reminiscent of a Ward 
identity, and is most probably a consequence of gauge symmetry, although we have not been 
able to demonstrate this convincingly. It can be checked that a corresponding relation holds 
between the coefficients \ and a in the RGE for W T [p] : 

2J dy± 5fr(y ± ) =a(X±) (427) 



(this is a rather obvious, consequence of eqs. ( 4.1S ) and ( 1.23a )) and, more significantly, also 



between the coefficients x an d cr of the original RGE (|1.16 ) for W r [p] 
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(4.28) 



This last relation is more significant in that it sheds some more light on the classical polar- 
ization term 5a : this is such as to make the two relations ( [4.27 ) and Q4.28|) consistent with 



each other, and with the p-dependence of the gauge rotation (|1.21a|) from x t° X- That is, 
5a has precisely the right value to compensate the functional derivatives of the Wilson lines 
V and V* in eq. flPIij) . 



There are at least two important consequences of the relation (|4.18|). The first one 



refers to the cancellation of infrared divergences in the evolution equations for gauge- 
invariant quantities. We shall see in Sect. |V B| that, when the Balitsky-Kovchegov equation 



is derived directly from the RGE in Hamiltonian form Q4.22j ) — where eq. Q4.18|) has been 



already used — , the IR divergences are automatically absent. But if one rather starts with 
the original RGE ( |4.1| ), where both rj and v are still present, then each of these terms will 
individually give rise to IR singularities, which will cancel only in the final sum. We expect 
similar cancellations to work for all gauge-invariant observables, but a general proof in this 
sense is still lacking. This would be the non-linear generalization of the standard cancellation 
of IR divergences between virtual and real corrections. 

The second consequence refers to the behaviour of the weight function W^a] at large 
"times" r, which is what governs the high energy limit of gluon correlations and, ultimately, 
of hadron cross-sections. To better appreciate the connection between eq. ( 4.1 8|) and the 



large time behaviour, it is useful to consider first the simpler example of the Brownian 
motion, which is governed too by a diffusion equation similar to eq. (|1.16| ), but for which an 
explicit solution can be readily worked out. Consider thus a particle suspended in a highly 
viscous liquid, and in the presence of some external force. The particle performs a random 
walk because of its random collisions off the molecules in the liquid. The relevant quantity 
the analog^ of W^ja] — is P(x, t) (= the probability density to find the particle at point 



x at time £), which obeys the following (Fokker-Planck) equation [|43"| , |4"4" 

where D is the diffusion constant (for simplicity, we assume this to be truly a constant, 
i.e., independent of x or t), and F l (x) is the "external force" (actually, the force divided by 
the viscosity). If F l = 0, the solution is immediate, and reads (with the initial condition 
p(x, o) = 5®(x)): 

p ^ = wk^ exp {~^}- (430) 

This is normalized to unity, as it should, at any t (J d 3 xP(x,t) = 1), and goes smoothly 
to zero at any x when t — ► oo (runaway solution). Thus, for sufficiently large time, the 



9 This analogy will be developed in more detail somewhere else [41]. 
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probability density is quasi- homogeneous (no x dependence), but it is everywhere close to 
zero (to cope with the normalization condition). The particle simply diffuses in all the 
available space. 

This situation changes, however, if the motion of the particle is biased by an external 
force. Assume this force can be derived from a potential: 

dV 

(4.31) 



dx % 



Then the time-independent distribution Pq(x) ~ exp[— /3V(x)] is a stationary solution to 
eq. ( |4.2iJ| ) provided (3D = 1. Of course, this solution is acceptable as a probability density 



only if it is normalizable, which puts some constraints on the form of the potential. But 
assuming this to be the case, then Po(x) ~ e _/3V/ represents an equilibrium distribution 
which is (asymptotically) reached by the system at large times. Once this is done, all the 
correlations become independent of time. 

Returning to our quantum evolution, a stationary distribution would correspond to 
a non-trivial fixed point of the RGE (|1.16|) (a solution Wo [a] which is normalizable and 
independent of time). If such a distribution existed, then the high energy limit of QCD 
scattering would be trivial: at sufficiently high energies, all the cross sections would be- 
come independent of the energy. The relation ( [4.18[ ) between the coefficients in the RGE 



guarantees, however, that this cannot happen: The effective force in eq. ( 4.1 7| ) vanishes, 
and the corresponding evolution Hamiltonian (|4.23|) is a purely kinetic operator, which de- 
scribes diffusion on the group manifold. We do not expect any non-trivial fixed pointQ for 
this Hamiltonian, just runaway solutions, and this is indeed the case for the approximate 
solutions to eq. flO^) found in Ref. (21). 



V. RECOVERING SOME KNOWN EQUATIONS 

If (0[a] ) T is any observable which can be computed as an average over a, as in eq. ([4.9D, 
then it satisfies an evolution equation given by ( [4.10|) , namely, 

j-(OH ) T = ( ~ . S . , 7fi . * N 0[a]\ , (5.1) 
Or' \2 5a«(x ± ) ' xy 5a b T (y ± ) / 



where we have also used the RGE in Hamiltonian form, eq. ( 4.22 ), and we have integrated 
twice by parts within the functional integral over a. In what follows we shall apply this 
equation for two choices of the operator O [a] : 

- The unintegrated gluon distribution function^ (i.e., the integrand of eq. Q1.4Q ): 



10 Of course, a constant distribution Wo = const, (no dependence on a) would be a stationary 
point for the Hamiltonian, but this should scale as 1/V with V = the volume of the manifold - 
to be normalizable — , and thus would vanish as V — > oo. 

n This is independent of k + since the electric field J- l+ {x) is almost a 5-function in x~; cf. eq. ([IT 
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ip T {kl) ee kl(\^(k\k ± )f) r (5.2) 

in the weak field approximation; we shall thus recover the BFKL equation, as expected. 
- The following two-point function of the Wilson lines (see below for its interpretation): 

S T (x ± ,y ± ) = (ti(V\x ± )V{y ± ))) , (5.3) 



for which we shall recover the Balitsky-Kovchegov equation (pifOJ . 



A. The BFKL limit 



The BFKL equation is obtained in the weak field (and source) approximation, in which 
one can expand the Wilson lines in the Hamiltonian to lowest non-trivial order. For instance, 
the BFKL Hamiltonian in the present formalism is obtained by replacing, in eqs. ( 4.23 ), 



1 - V}V X — ► ig(a a (x ± ) -a a (z ± ))T a , a a (x ± ) ee J dx~ a a (x~, x±) . (5.4) 



Clearly, this is formally the same as the perturbative expansion (i.e., the expansion in powers 
of g) of the Hamiltonian to lowest order. After this expansion, the Hamiltonian takes the 
generic form 

Hbfkl ~ «r a 7"' ( 5 - 5 ) 
oa da 



which is like the second-quantized Hamiltonian for a non-relativistic many-body problem: 
In this weak field limit, the quantum evolution is diagonal in the number of fields, i.e., it 
couples only correlators (a(l)a(2) • • -a(n)) T with the same number n of fields. The BFKL 
equation is the corresponding evolution equation for the 2-point function (|5.2| ). 

To write down this equation, we first observe that, for such weak fields, the solution 
to the Yang-Mills equations ( |1.3| ) can be linearized as well: 

F?(k) « HVll&pJJz), (5.6) 



which implies an approximate expression for the unintegrated gluon distribution (|5.2j ): 

<Pr(kl) « ( Pa (k ± )p a (-k ± )) T , (5.7) 

where 

p a (k ± ) = J d 2 x±e~ ik± ' x± p a (x±), p a (x±) ee J dx~ p a (x~ , x±) . (5.8) 

We thus need the equation satisfied by the charge-charge correlator (pp) T , which follows 
quite generally from the RGE (|1.16|) in the p-representation: 
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-^;(pa(x±)pb(y±)) T = a s (o- a (x±)pb(y±) + pa(x±)cr b (y±) + Xab(x±,y±)) T - (5.9) 



For generic, strong, fields and sources, the r.h.s. of this equation involves n-point correlators 
(p(l)p(2) • • • p(n)) T of arbitrarily high order n. But in the weak field limit, where a is linear 
in p and x is quadratic, this becomes a closed equation for the 2-point function, in agreement 
with the general argument after eq. ( |5.5| ). The coefficients x an d cr in this limit will be now 
obtained by expanding the general formulae ( p. 26 ) and ( 3.1 9|) to lowest order in p. 



Consider x first. To the order of interest in p, one can replace, in eq. fl2.26| ), 



AUx±) « -d*a(x ± ) « (dyV 2 ± )p(x ± ), « d\ (5.10) 

which after simple algebra yields (in matrix notations: (p x p y ) ab = pTp'vi e ^ c : 



Xab( X ^V±) 



'x ry 

d 2 p ± e ip ±< x ±- y± ) 



(2vr) 2 p\ 

X {pxPy - ip x 



or, after a Fourier transform, 

Consider similarly : to linear order in p, this can be extracted from any of the expressions 
Q3.19] ), ( |3.21| ), or ( |3.28| ), so let us choose the latter expression, for convenience. To this aim, 
it is enough to replace, in eq. ( |3.32|) , 



Tr(T a VX) » tgN c [a a (x ± ) - a a (z ± )}, (5.13) 
as appropriate for weak fields, and thus get, after simple algebra, 

°^ = - N ^j^m£^?- (5 - 14) 



By inserting eqs. Q5.12Q and (|5.14| ) into the evolution equation (|5~^ ), and using ( |5.7D , 
one finally obtains: 



dip T {k\) a s N c 



f k 2 
J vAk\-v\Y 



,KX ' r d 2 Pl _ 2n kl ,^ T (ki), (5.15) 



dr it 2 J p 2 ±(k± — p±) 

a s N c r 2 k\ 
2n 2 J P± p 2 ± (k ± -p ± ) 2 

which coincides, as anticipated, with the BFKL equation [O. The first term in the r.h.s., 
which here is generated by x > is the real BFKL kernel, while the second term, coming 
from o"(°) , is the corresponding virtual kernel. 
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B. The Balitsky Kovchegov equation 



Note first that, although written in a specific gauge — namely, the co variant gauge 
where the color field of the hadron reads = 5^ + a a {x) — , the quantity ( |5.3|) has a gauge 
invariant meaning, as the COV-gauge expression of the gauge-invariant operator 

S T (x±,y±) = (trL(x±,y ± ))^ (5.16) 

where L(x±,y±) is the Wilson loop obtained by closing the two Wilson lines along x~ (at 
x± and y±, respectively) with arbitrary paths joining x± and y± in the transverse planes at 
x~ = oo and x~ = — oo. 

Physically, S T (x±, y±) is the S-matrix element for the scattering of a "color dipole" off 
the hadron in the eikonal approximation. To clarify this, consider DIS in a special frame 
(the " dipole frame" ) where most of the energy is put in the hadron which moves along the 
positive z direction — in this respect, this is like the infinite momentum frame, so one can 
use the CGC effective theory for the hadron wavefunction — , but the virtual photon is itself 
energetic enough for the scattering to proceed as follows: The photon first fluctuates into 
an energetic quark-antiquark pair (a "color dipole" ) which then propagates in the negative 
z (or positive x~) direction, i.e., towards the hadron, against which it scatters eventually, 



with a cross section (see, e.g., |I5|,20|,ID[ for more details) 



cr dtp oie(r,r ± ) = 2 J d 2 b ± ^(tr(l - V^(x ± )V(y ± ))) ^ (5.17) 

Here x± and y± are the transverse coordinates of the quark and antiquark in the pair, r± = 
x±—y± is the size of the dipole, and b± = (x±+y±)/2 is the impact parameter. VJ and V y are 
Wilson lines in the x~ direction, defined as in eq. ( |1.8|) but in the fundamental representation. 
They describe the eikonal coupling of the quark and antiquark to the color field = a a 
of the hadron. The average over the hadronic target in eq. ( |5.17| ) is here understood as 
an average over a in the sense of eq. ( |4.9|) . When the scattering energy increases (i.e., r 
increases), the weight function W r [a] for this average changes due to quantum corrections, 
and so does the scattering cross section. This change is governed by eq. ( |5.1|) with O[o\ = 
S T (x±,y±), that we shall construct now explicitly. 

The final outcome of this calculation, to be detailed in Appendix D, is an evolution 
equation for S T (x±,y±) that we write down here for a color dipole in some arbitrary repre- 
sentation R, since this is not more difficult: 

§p(*xMV v )) r = -J / ^ {x J^ y f_ z±) 2 (Cr*rMV v ) - tr R (V^V z V^V y )) T . 

In this equation, all the Wilson lines V\ V, and the color group generators t a , t b belong to 
the generic representation R, and Cr = t a t a . 

In particular, for the fundamental representation, the following Fierz identity: 



(5.1i 
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t iUkl = Igilfijk '- 1 § ij § kl / 5 19 ) 

a a 2 2N C y ' 

allows one to simplify the color trace in the last term in eq. (|5.18|) . By also using Cr = 
(iV| — 1)/2N C , one finally obtains (with R = F kept implicit) 

^HvX)) T = / d2 ^ (x J^ y f_ z±f (NMvX) - tr(Kfv.)tr(vJy tf )) T , (5.20) 

which coincides with the equation derived by Balitsky within a rather different formalism 



18| : via the operator expansion of high-energy scattering in the target rest frame. 



The first observation is that the above equations are not closed. They relate a 2- 
Wilson-line correlation function to a 4-line function, for which one can derive an evolution 



equation too [T8|], but one may already guess that this would be not the end of the story: 
The 4-line function will be in turn coupled, via its evolution equation, to a 6-line function, 
and so on, so that one is really dealing with an infinite hierarchy of coupled equations of 



which eq. ( |5.20| ) is just the first equation. 



But a closed equation can still be obtained in the large N c limit, since in this limit the 
4-line correlation function in eq. ( |5.20| ) factorizes: 

tr(y2V,)tr(V2V y j) — ►(trCK+V,)) (tr(V z %)) for N c - oo, (5.21) 



and eq. (|5.20|) reduces to a closed equation^ for the quantity Af(x±, y±) = (tr(l — V£V y )) T 
(= the forward scattering amplitude of the color dipole with the hadron) : 

^Ky = ^J ^ {x± }^ y f_ z±) 2 + Ky ~ Ky) ~ K z M Z y) ■ (5.22) 

This coincides with the evolution equation obtained by Kovchegov within still a different 
formalismF], namely the dipole model by Mueller [26 . 



The second observation is that the Wilson lines in the adjoint representation which 
were originally present in the Hamiltonian (|4.23|) have melted together with the Wilson lines 
in the generic representation of the external dipole to generate the color traces in the r.h.s. 
of eq. ( |5.18| ). This melting has been made possible by relations like: 



12 



A different closed equation has been recently obtained in Ref. p8| by evaluating the 4-line 



correlator in eq. (5.20) with a Gaussian weight function for finite N c . Note however that the 
general solution to the RGE ( |4.22| ) is not a Gaussian |41] and, moreover, eq. ( |5.20| ) is not even 
consistent with a Gaussian, or "mean field", approximation to the general RGE pllfji 



13 A similar non-linear equation describing the multiplication of pomerons was suggested in Ref. 
[j|] and proved in P] in the double-log approximation. More recently, Braun has reobtained this 
equation by resumming "fan" diagrams 
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V} b t h = Vt a V\ (5.23) 

which allows one to trade a group element in the adjoint representation (that we denote 
here with a tilde, for more clarity) for a pair of group elements V and in the generic 
representation with generators t a , t b . Note also the simple correspondence between the two 
terms within the braces in the r.h.s. of eq. (|5.18|) and the four terms in eq. (|4.2|) for rj : The 



first two terms, 1 and VJV y , in eq. ( |4.2j ) have contributed equally to the first term, involving 
(trR(VjiVy)), in eq. fl5.18| ), while the other two terms, VJV Z and V] V y , have given identical 



contribution to the remaining piece, involving the 4-line function, in eq. (|5.18|) . 

The third important observation refers to the transverse kernel in eq. (|5.18|) , or G5.20|) . 
A brief comparison shows that this is not the same as the original kernel )C(x±, y±, z±), 
eq. Q4.16D , of the RGE. Rather, this has been generated as (see the explicit calculation in 
Appendix C) 

(x_\_ — y±)^ 

K.(x_ L ,x±,z±) +)C(yx,y±,z±) - 2JC(x±,y±,z±) = ^ . (5.24) 

{x± - z±) 2 {y± - z ± ) 2 

The final result, also known as "the dipole kernel" (since it is proportional to the probability 
for the decay of a dipole of size x± — y± into two dipoles of sizes x± — z± and y± — z±, 
respectively), has the remarkable feature to show a better IR behaviour than the original 
kernel (|4.16| ). When z± ^> x±, y±, the r.h.s. of eq. Q5.24Q decreases like (x± —y±) 2 /z±, so its 



integral / d 2 z± is actually finite. This guarantees, in particular, that the evolution equation 
( 5.22j ) in the large N c limit is free of IR problems, and it strongly suggest a similar property 



for the exact equation ( 5.2U| ), although a rigorous proof in this latter case would still require 
a knowledge of the large-z_i_ behaviour of the 4-line function (tr(V^V^)tr(V^V^)). (All of 
the kernels occuring in the evolution equation for this 4-line correlation function are of the 
dipole type too flIE| .) 

The dipole kernel (|5.24j ) has ultraviolet poles at z± = x± and z± = y±, but these are 



inocuous, since the accompanying Wilson-line operators in eqs. ( |5.18| ) or ( |5.20| ) cancel each 
other at these points. 



VI. CONCLUSIONS AND PERSPECTIVES 



In this paper we have completed the construction of the renormalization group equa- 
tion for the Color Glass Condenstate by explicitly computing the coefficients in this equation, 
the one-loop real and virtual contributions, to all orders in the color glass field. The resulting 
expressions appear to be quite simple — given especially the complexity of the calculation - 
and to exhibit some remarkable properties with important consequences. We have started 
the exploration of these properties and of some of their consequences, but there is much room 
left for further studies. In particular, the deep reason for some of these properties, like the 
functional relation between the real and virtual contributions, is not yet fully understood. 

Our results turn out to be equivalent, at least as far as the evolution of Wilson line 
operators is concerned, with previous results obtained by Balitsky |1| and (in the large 
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N c limit) by Kovchegov within quite different approaches: the operator expansion of 



high-energy scattering amplitudes in Refs. [ lql , and the Mueller dipole model in the case of 



Ref. [ 19 1 . In fact, our RGE ( f4.1| ) appears to be the same (up to some change of variables) 
as the functional equation written down by Weigert to summarize in closed form Balitsky's 
infinite hierarchy of equations This agreement between our respective results — in spite 
of the profound differences between the various formalisms staying behind them — gives 
us even more confidence that these are indeed the right equations to describe non-linear 
evolution at small x in the leading logarithmic approximation. 

On the other hand, our results appear to disagree from those reported in previous 



attempts [|17| to compute the coefficients in the RGE. There are several noticeable differ- 
ences between our calculations and those in Refs. |L7| - chiefly among them, the different 
constructions of the background field propagator (especially in relation with the gauge-fixing 
prescription), the different ways to treat the longitudinal extent of the color source, and the 
role of the classical polarization term, which has been missed in Refs. |T7[ - which may 



explain the different results, but a deep understanding of these differences requires more 
study. 

Now that the RGE is known explicitly, the next objective is, of course, to try and 
solve it. Significant progress has been already achieved in this direction, with interesting 



consequences. In Ref. ||21|| , some approximate solutions have been obtained in a mean field 
approximation in which the coefficients in the RGE have been replaced by their expectation 
values, self-consistently determined from the ensuing approximation to the weight function. 
The results thus obtained distinguish between two regimes: 

- A short-distance, or high (transverse) momentum regime k± ^> 1/Q s (r), where one 
recovers the original Ansatz of the McLerran-Venugopalan model, namely that of a system 
of independent color sources whose density increases with r. This is a low density regime, 
where the linear evolution equation is valid. 

- A long-distance, or low momentum (k± <C 1/Q s (t)), regime, where the system of color 
charges rather forms a Coulomb gas with long-range correlations in the transverse plane. 
This is a high-density regime which exhibits gluon saturation, with a maximal gluon phase- 
space density of order l/a s . 

The separation scale Q s (r) between these two regimes (the "saturation scale") is esti- 
mated as |IJ,0,^,^]]] (with R = the hadron radius) 

R 2 Q 2 s (r) ~ a s j^-j xG(x, Q 2 s (r)) oc e Ca - , (6.1) 

where the last estimate (with C a pure number) is based on the linear evolution equation in 
the double logarithmic approximation and with a fixed coupling constant a s . If one rather 
uses the one-loop running coupling constant a s (Q 2 ) with Q 2 = Q 2 (t), then the saturation 
scale rises more slowly, as Q 2 s {t) ~ e^~, with b some number. 

It is interesting to note that these predictions of the RGE are in a good agreement 
with recent investigations of the Kovchegov equation ( [5.22| ) , for which various solutions have 



been obtained via a combination of analytical and numerical methods [|l9| ,p^-|24]| . Although 



the large N c limit (in which Kovchegov' equation becomes applicable) is a priori not the 
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same as the mean field approximation of Ref. f2~I| , it is nevertheless clear that these two 



approximations are closely related, which explains the similar results. Thus, more effort is 
required to relax these approximations, and compute 1/N C corrections. A possibility to do 
so is via lattice simulations. The exact solution to the RGE (|4.22|) can be given a path 



integral representation, to be presented somewhere else [ |1 1 . Another path integral has been 
recently proposed by Balitsky [[4^], and provides an exact solution to his coupled equations. 
All these path integrals are those of a quantum field theory in three Euclidean dimensions 
[two transverse coordinates plus "imaginary time" = r (rapidity)], and thus are well suited 
for lattice calculations. 

A better control of the solutions to the RGE would allow one to extract its phe- 
nomenological consequences, notably in relation with high-energy deep inelastic scattering 
and ultrarelativistic heavy ion collisions. Possible signatures of saturation in DIS and heavy 
ion collisions have been already discussed in the literature p9|-[38|j9|]. It would be interesting 



to see how these and related predictions are affected by the details of the effective theory at 
small x. 
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APPENDIX A: THE BACKGROUND FIELD GLUON PROPAGATOR 



The propagator G^(x,y) of the semi-fast gluons in the background of the tree- level 
field ( |1.5| ) and in the LC-gauge (a+ = A£ = 0) has been obtained in Sect. 6 of Paper I (see 
also Refs. |39| , ^0| ) via the following three steps: i) First one has constructed the background 



field propagator G(x,y) of a scalar field, ii) Then, one has computed the gluon propagator 
G^(x, y) for a quantum gluon in the temporal gauge a~ = (but with the background field 
= b^ % A % c still in the LC-gauge). Hi) Finally, the LC-gauge propagator has been 
obtained from G^ b via a gauge rotation. 

Below, we shall present the final results of these calculations. For convenience, we 
start with the free propagator. 



1. The free LC-gauge propagator 



Consider first the free propagator (no background field) GQ V ab = S a bGQ P '. With the 
retarded prescription advocated in Sect. 1.3, this reads, in momentum space: 
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Gi(p) = G Q (p), G?(p) = ^— G (p), 

G#(p) = ^Go(p), C --(p) = PV |f Gob) , (Al) 

where Go(p) = l/(2p + p~ — pi + ie) [the free propagator of a massless scalar field], and PV 
denotes the principal value prescription: 

pv'.if ' + '). (A2) 
p+ 2 \p+ -ie p + +iej 

Note that, strictly speaking, the gauge-fixing prescription in eq. (|A1|) is "retarded" only as 
far as the component G ~ is concerned; by hermiticity, the corresponding prescription in 
Gq 1 is advanced, while in Gq ~ it is a principal value. 

Other prescriptions that are referred to in the main text are the "advanced" prescrip- 
tion (see, e.g., Refs. [p,14,16|) which is obtained by changing the sign of ie for the axial poles 
in eq. ( |A1| ) [that is, this is advanced for the component G ~], and the PV-prescription, for 
which l/p + = PV(l/p + ) in all the components of the propagator. 

For further reference, we note also the expression of the scalar propagator in the x~ 
representation: 

G (x-,p-, P± ) = e-**-*- n _^__ — 2 



2p 



2n 2p + p~ — p ± + ie 

{0(aT)0(p-) - 6(-x-)6(-p-)} e l ~ x . (A3) 



2. The background field propagator of a scalar field 



Since the background fields are static, and thus we have homogeneity in time, it is 
convenient to work in the ^"-representation, that is, to construct the propagator G(x, y,p~) 
for a given p~. This also makes it easy to implement the strip restriction ( p..31| ). The 
propagator is known only for the discontinuous background field A l (x) = 6{x~)-Vd l V^ in 
eq. flOl), which is how the actual field of eq. ( |1.5| ) is "seen" by the semi-fast gluons^- The 
final result for G(x, y,p~) can be written as: 

G(x,y,p-) = G (x - y,p-){9(x-)9(y-)V(x ± )V\y ± ) + 6(-x-)6(-y-)} 
+ 2ip~ J d 3 zG (x- z,p~) 8{z~) G Q (z — y,p~) 

x {6(x-)6(-y-)V(x ± )V^z ± ) - 6(-x-)6(y-)V(z ± )V\y ± )} . (A4) 



14 With p~ constrained as in eq. ( |1.31 ), the typical momentum scale p + ~ P±/p~ for longitudinal 
dynamics is relatively small, p + <C A + , so that the semi-fast gluons are unable to discriminate the 
internal structure of the -function in eq. dl 
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It can be easily verified that this function is continuous at both x =0 and y =0. 

On eq. (|A4j) , we distinguish two type of contributions, "crossing" and "noncrossing" , 
corresponding to trajectories which cross, or do not cross, the plane at x~ = where 
is located the singularity of the background field (or its source). Consider "noncrossing" 
trajectories first, for which x~ and y~ are of the same sign: When they are both negative, 
the propagation takes place in a domain where the field vanishes; this is therefore free 
propagation. When x~ and y~ are both positive, the gluon propagates in a background field 
which is just a gauge rotation, A % = -Vd % V^ ; thus, the net effect comes from the difference 
between the gauge rotations at the end points. If the trajectory crosses the discontinuity 
at x~ = 0, say in going from y~ < to x~ > 0, then there is also a gauge factor V\z±) 
associated with this crossing (at some arbitrary z±). 

We thus write, with obvious notations, G = G^ 1 ' + G^ c '. In the crossing piece G^ c \ it 
is further possible to replace: 

9{x-)6(-y-) — 0(p-), e(-x-)9(y-) — 9(-p-), (A5) 

because of the correlation between the sign of p~ and the direction of propagation in x~, as 
manifest on eq. flA3|) . 

3. The gluon propagator in the temporal gauge 

This has the following non-zero components (G = the scalar propagator in eq. (|A4|) ) : 

G^(x,y,p-)=5^G(x,y,p~), 

G +i (x,y,p-) = —ViG(x,y,p-), G l+ (x,y,p-) = - ± G(x, y,p~) 2?Ji , 
p p 

G ++ {x,y,p-) = jl- {v l x G(x,y,p-)Vl l + 5®(x -jfi} . (A6) 



Note that, because of the strip restriction ( |1.31|) on p , the operator 1/p is never singular. 

In the above equations, T> % = d l — igA 1 and = + igA-* with the derivative 
acting on the function on its left. By using the Wilson lines in eq. ( |A4|) , it is possible to 
convert these covariant derivatives into ordinary derivatives. This relies on the following 
identities, valid for any function 0(x) (cf. eq. (|1.5[)) : 



V 



U{x)0{x)] = U{x)d l O{x), [0(x)U\x)]v ]l = (9 i O(x))?7 + (f). 



(A7) 



This gives the following expressions for the crossing and non-crossing pieces of, e.g., G ++ : 

2i f 

-- — J d 2 z± d l x G (x~,x ± - z±) dlG (-y~,z ± - y±) 
x {e(p-)V(x ± )V\z ± ) - 6(-p-)V(z ± )V\y ± )} , (A8) 



C++( C ) ( 



x,y;p~ 
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and 



G ++in \x,y;p-) = {6{x-)6{y-)V(x ± )Vi{y ± ) + e{-x-)e(-y-)} 



x^^a*G (f-j/;p-) + t {3) (x-y)j^y- (A9) 

In the crossing piece we have also performed the replacement ( |A5|) . 

In the zero field limit, the above equations yield the free temporal-gauge propagator 
in the expected form: 

G l o J (p) = S ij G Q (p), Gi+(p) = Gp{p) = £ G (p), G+ + (p) = ^ G (p) , (A10) 
vwhere G (p) = l/(2p + p~ — p\ + ie). 



4. The gluon propagator in the LC gauge 

This is related to the temporal-gauge propagator presented above via the relation: 
iG£{x,y) = (Taf^d,)) = (T(d» -V^d+)l(a» -V»±d + Y) , (All) 



or, more explicitly, 



G ij = qH G +j + G i+ —T)ti - G ++ ^-V^ , (A12a) 

G - i = _|_G + ._|_G++±pt. i (A12b) 

(A12c) 

G-=- 8 -G++ 8 -, (A12d) 

where we use the convention that the derivatives written on the right act on the functions 
on their left; e.g., d~F d~ = d~d~F(x, y). 

The axial poles at p + = in the expressions above are regularized according to the 
"retarded" prescription in eq. (|A1|) . This implies that the operators l/d + written on the 
left (right) of G^ u in eq. (|A12|) are regularized with a retarded (advanced) prescription. For 
instance, 



1^1 1_ ^ 1 



where 
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1 



[X 



\y-) 



dp + e ip+( - x - y ) 
2n p + + ie 
dp + 



In p + — ie 



-iv [x — y 



id{y- 



(Al4a) 
(A14b) 



In going from eq. ( All ) to eq. ( |A12| ), we have also used : 

(x 



(A15) 



APPENDIX B: ALTERNATIVE CALCULATION OF x 3 

In this Appendix we shall verify that a complete calculation of %3 using the full tem- 
poral gauge propagator G ++ leads to the same result Q2.13Q as the simplified calculation 
with G ++ replaced by its free counterpart Gq + . 

According to its definition in eq. (|2.1| ), G ++ involves the following matrix element (the 
equal-time limit x + = y + is understood) 

M(x ± ,y ± ) = (x- = 0,x ± \-^G ++ -^\y- = 0,y ± ) = + (Bl) 

io R %a A 

with both crossing (Ai^) and non-crossing (Ai^) contributions. The crossing piece is 
evaluated with the propagator G ++< - c - ) in eq. (|A8|), which yields 

J strip 2tt J Z7T J 2tt p + + ie k + — ie 

dp~ 2% f dp + f dk + 1 1 f , , . 

^ — - ^~ ^ — i~r~ it — ~ / dT± ^ ' k± ) 

strip 2tt p J 2tt J 2tt p + + ie k + — ie J 
G (p)G (k) {6(p-)V(x ± )V^(z ± ) - 9(-p-)V(z ± )V^y ± )} . (B2) 

(We have also used here the notation (|3.9|) for the integrals over the transverse phase-space.) 
This vanishes since the integrals over p + and k + generate mutually excluding 9 functions : 

/fr^ G °M^-^< /fr^^w^^- (B3) 

Thus, M (c) = 0, in agreement with the arguments following eq. 

As for the non-crossing contribution, this is most easily evaluated as (cf. eqs. ( |2.4| 
and QA§D) 

M {n) = f dz- x [ dz 2 e(-z;)6(-z 2 )G ++ ^(z;,x ± ,z-, y± ) 



46 



x / e'^-<^-»J (p'i G„(p) + l) 



(2tt)* 



G() (p) (B4) 



The integrations over z l 2 are performed using 



dz~9(-z~)e Tip+z ~ = ±l (B5) 

p + =F ie 



which gives 



M^= \%- \ d 4- [P^ e^(—> -1- ^ -1- G (p). (B6) 
7 2vr i 2vr i (2tt) 2 p+ + ie p- p+ - ie V ; 

This is the same expression ( |2.8|) as obtained by using directly the free propagator Gq + . 
From now on, the calculation proceeds as in Sect. [II B| and leads to eq. ( p. 13 ) for X3, as 
anticipated. 



APPENDIX C: MORE DETAILS ON THE VIRTUAL CORRECTION 

In this Appendix we collect some calculations pertinent to the virtual correction a(x) 
which have not been included in the main text. 

(a) Computing &2 

We first compute the second contribution <7 2 to eq. (|1.36|) for the induced source and 
show that this is indeed a pure tadpole, as stated in Sect. 1 LI A|. This is the quantum 



expectation value of the second piece in eq. ( |1.27|) for 5p^\ that is: 



2 

= -*jj-J d y + J dz + M abcd (x + ,y + ,z + ) Pb (x)G~ cd -(y + ,x;z + ,x), 

M abcd (x + , y + , z + ) = 9{x + - z + )9(z + - y + )Tr(T a T b T c T d ) + d(z + - y + )d{y + - x + )Tr (T a T c T d T b ) 
+ 9(z + - x + )6(x + - y + )Tr{T a T d T b T c ). (CI) 

By using eq. (|A12d|) for G together with time homogeneity, one can write: 

G«(y + ,2;z + ,Z) = J ^- e ~"- '"*-'*> G;i(x,x;p-j 

t^^^^lk^lk^ <C2) 
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In eq. (|CID , this is multiplied by Pb(x), so one can replace x~ ~ 0. Then ( |C2"D is essentially 
the same matrix element as in eq. (|2.4j ), except that this is now evaluated for arbitrary 
time variables y + and z + . But the arguments in Sect. [LI Aj, according to which this matrix 



element can be equivalently evaluated with the free propagator G 0cd = 5 c dG Q , are still 
valid, since they rely just on the longitudinal structure of the propagator together with our 



axial prescriptions. Thus, for a; =0 eq. (|C2j ) is the same as (cf. eqs. (|2.8|) - (|2.10|) ) 



5 cd G^(y + ,0-,x ± ;z + ,0-,x ± ) = -i5 cd [ ^-e^^-^ p~e(p-) / 

J strip Z7T J 



d 2 P± 1 



(2tt) 2 - 2 ' 



(C3) 



which is the same matrix element that would enter at BFKL level. Because of that, also 



the final result for &2 is the same as at BFKL level, and is therefore well known [13]. For 
completeness, let us nevertheless finish its calculation. 



With the trivial colour structure of eq. (|C3j ), the colour traces in eq. (|C1|) become 
st raight forward : 

M abcc {x + ,y + , z + ) = [e(x + - z + )6(z + - y + ) + 6(z + - y + )6(y + - x+)}Tr(T a T b T c T c ) + 

+ 6(z + - x + )6{x + - y + )Ti{T a T c T h T c ) 
= 6{z + - x + )6(x + - y + )Ti{T a [T c ,T b ]T c ) + 6(z + - y + )Tr{T a T b T c T c ) 

y+). (C4) 



-S" b -SLe(z + - x + )6(x + - y + ) + 5 ab N 2 c e(z^ 



In going from the first to the second line above, we have used the identity: 

6(x + - z + )6(z + - y + ) + 6(z + - y + )6(y + - x + ) = -6{z + - x + )6{x + - y + ) + 0(z + - y + ) . (C5) 



To obtain <5"2, one still has to perform the integrations over y + , z + and p~ . After these 
operations, we get a non- vanishing contribution only from the first term (the term involving 
two ^-functions) in the last line of eq. (|C4|) . This involves: 



dy + / dz + 6(z + - x + )6(x + - y + ) e 



-ip (y+-z+) 



-1 



(p- 



(C6) 



where the prescription at p~ = is irrelevant since p~ is anyway restricted to the strip ( |1.31| ). 
This restricted integration over p~ then generates the expected logarithmic enhancement, 
cf. eq. (|2.11| ), so that the final result reads: 



a 2 {x) 



9 2 N C 
2tt 



\n l -p a (x) J 



d 2 p ± 1 
(27ppi' 



(C7) 



which is linear in p, and a pure tadpole, as anticipated. 

(b) The "non-crossing" contribution to ai vanishes 

Let us now check that the "non-crossing" contribution to eq. ( |3.7| ) for &i vanishes 
indeed. This involves the propagator G ++( - n ^ of eq. (|A9|) , and is therefore proportional to 
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a$ oc VI f dz~ 6(x- - z~) G ++ ^(z',x-;x ±J y ± ) Pj* (C8) 
' dz~ 6{x~ - z-)Vi{0{z-)9{x-)V x V^ + 6{-z-)9{-x-)}(didiG ) V%\ x= 
dz~ 6(x- - Z -){9(z-)0(x-) + +6(-z-)6(-x-)}v{Go(z- - x~,x± - y ± ) 



x±=y± 



which is trivial in color (oc 5 ab ), and therefore gives no contribution after taking the color 
trace with T a (cf. eq. ( |3.2| )). In going from the second to the third line in the equation 
above, we have used (cf. eq. (JO)) 



VI = 9(-x-)dl + e{x-)Vl, V i oa = d i - igAU^) = & + V&V\ (C9) 
to deduce that, for any function 0(x±), 

v^iyo) = v&o, (ov j )v^ = (d l o)v ] . (cio) 

This has allowed us to convert the covariant derivatives into ordinary derivatives acting 
on G • Finally, in the local limit x± — > y±, V X V^ — > 1, and all the color structure has 
disappeared. 

APPENDIX D: DERIVING THE BALITSKY KOVCHEGOV EQUATION 



We present here some of the calculations leading to eq. (|5.20 ) in the main text. The 



starting point is eq. ( |5.1| ) with 0[a] = S T (x±,y±) given by eq. (|5.3|). This gives (the colour 



group representation R is arbitrary): 

d I. \ 1 r ,o r ,o / 5 „ h 5 



where (cf. eqs. (gj ) and (|4.16|) ) : 



in' 



The first derivative in eq. ( |D1|) is immediately performed with the help of eq. ( 4.12|) (with 
tr = tr R ) 

tr(t£V w ) = z<?tr(*VJ^)(5 OT -<W, (D3) 



5a b T (v ± ) 



where 5 XV = — i>j_). To perform similarly the second derivative in ( pi|) , we consider 

separately the contributions of the four terms within rj^, eq. ( [D 2| ) . 

The first term, involving the unit matrix 5 ab , is immediate: 
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■^x-S7TT tr (Kt^) = -9 2 Ctr(V^V y )(5 xv -8 yv )(5 xu -5 yu ). (D4) 



We have used here t a t a = Cr together with the cyclic invariance of the trace. 

For the second term, involving (V^V V ) , we use the identity to (|5.23|) to first write: 

(V^V v ) ab tr(t b VX) = (Vtrtr(V v k c V v VjV y ) = (V^rtr(V^V y ), (D5) 

where the second identity holds in the presence of the factor (6 XV — S yv ). The first factor 
{V^) ac in this equation commutes with the derivative 5/5a®(u±), since (T a ) a b = (this is 
similar to eq. (|4.19|) ). Thus, the contribution of the second term in eq. ( P2|) can be evaluated 



as: 



ig ($» - 5 yv ) (V^r jj^ tr(V^V y ) = -g 2 (5 XV - 8 yv ) (5 XU - 8 yu ) (V^ ac tr(t a V^V y ) 

= -g 2 C R tv(V^V y ) (8 XV - 8 yv ) (8 XU - 5 yu ) , (D6) 

where in writing the second line we have used eq. ( [5 . 2 3| ) once again, together with the 5- 
functions at u = x and u = y. Remarkably, the final result above is exactly the same as the 
first contribution in eq. (171). 



Consider similarly the third term (— V^V z ) ab in eq. ( |D2| ). It is easy to check that this 
commute with the derivative S/Sa^(u±) (as in eq. ([4.19|) ), so its contribution reads: 

d ^u^tJt—MVX) = -ig (8 XV - 5 yv ) (VtV z y b -£— tr(t b VX) 



9 2 (S xv - 5 yv )(V^) ac V z cb [tr(t b t a V}V y )6 xu - tv{t a t b V$V y )5 % 
g 2 (S xv - S yv ) (5 XU - 6 yu ) V z cb t<t b V^t c V y ) , (D7) 



after by now standard manipulations. 

A similar calculation shows that last term {—V}V v ) ab in eq. ( |D2| ) gives exactly the same 
contribution as in eq. ( |D7| ). By putting together the previous results, it is now straightfor- 
ward to obtain eq. ( |5.18p . In particular, the two integrals over u± and v± in eq. ( pip are 
trivially performed with the help of the ^-functions (S xv — S yv )(5 xu — S yu ), thus generating 
the dipole kernel ( |5.24|) . 
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